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1 .  INTRODUCTION 


The  Air  Force  Geophysics  Laboratory,  Optical  Physics  Division 
(AFGL, OPT.)  has  been  developing  a  high  resolution  radiative  transfer 
code  for  direct  point-to-point  transmittance  and  radiance  in  the 
earth's  atmosphere.  This  line-by-line  spectral  code  uses  the  AFGL 
line  atlas  and  an  algorithm  called  HIRACC^  to  construct  accur¬ 
ate  radiative  absorption  properties.  To  permit  computationally  fast 
manipulation  of  the  Voigt  lineshape,  an  atmospheric  transmittance  and 
radiance  model  was  developed  called  FASCODE"\  A  version  of  the  resulting 
code  along  with  certain  external  geometry  and  model  atmosphere  components 
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was  released  for  public  use  and  called  FASC0D1B  .  This  report  develops 
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and  documents  (incorporating  work  presented  in  Annual  Reports  1  and  2  *  ) 
modules  and  subroutines  for  the  generic  code  called  FASCOD  which  will  event¬ 
ually  be  merged  by  AFGL  into  a  integrated  code  to  be  designated 
FASC0D2 . 

The  AFGL  low  resolution  (typical  resolution  of  about  20 
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wavenumbers)  radiative  transfer  code,  LOWTRAN  5  ,  contains  extensive 
and  well  documented  gaseous  and  aerosol  model  atmospheres.  For  con¬ 
sistency  and  certain  possible  intercomparisons,  these  same  model 
i..tmospheres  are  built  into  FASCOD  and  this  process  is  discussed  in 
Seel , , n  2.  The  spherical  geometry  with  refractive  bending  at  the 
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layer  bovndaries  that  has  been  used  in  the  past  by  AFGL  in  LOWTRAN  5  , 
h«H  now  been  modified^  by  AFGL  to  include  continuous  refractive 
bendlin'  and  a  more  accurate  optical  path  routine.  As  part  of  Section  2 
this  new  geometry  package  along  with  zeroth-order  aerosol  scattering 
and  absorption  are  included  in  FASCOD.  The  modifications  presented 
in  Section  2  used  an  AFGL  version  of  FASCOD  called  FASCOD1A, 
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and  the  resulting  intermediate  code  is  called  FASC0D2-80  in  this  report. 

For  altitudes  above  40  to  60  Km,  certain  molecular  species 
can  be  found  to  have  nonequilibrium  population  distributions  in  their 
vibrational  states  due  to  photo-driven  chemistry  and  a  reduced  collision 
rate.  Consequently  the  usual  assumption  of  local  thermodynamic  equilib¬ 
rium  (LTE)  that  leads  to  the  Planck  function  and  absorption  coefficient 
used  for  thermal  radiation  sources  is  not  valid.  Degges  and  Smith^ 
have  modeled  such  phenomena  for  selected  vibrational  transitions  and 
have  predicted  non-LTE  population  densities  for  certain  atmospheric 
conditions  like  daylight  or  night  by  season  and  location  in  the  atmos¬ 
phere  (sun  angle).  Section  3  develops  the  necessary  theory  to  cal¬ 
culate  non-LTE  transmittance  and  radiance  and  modifies  FASCOD  to  ex¬ 
ecute  in  this  mode  when  non-LTE  population  densities  are  supplied. 

Code  changes  are  outlined  and  sample  computational  results  presented. 

Near  band  heads  where  the  far  wings  of  many  superposed  lines 
form  an  absorption  continuum,  the  Lorentzlan  line  shape  does  not  pre¬ 
dict  correct  results.  This  problem  of  the  sub-Lorentzian  line  shape 
is  discussed  in  Section  4  in  the  context  of  past  work  on  th  oretical 
modeling  and  experimental  data.  A  new  sub-Lorentzlan  form  factor  is 
developed  that  is  compatible  with  the  algorithm  for  line  modeling  in  FASCOD  and 
produces  continuum  absorption  predictions  for  CO2  that  are  as  accu¬ 
rate  as  those  using  more  complex  form  factors.  The  resulting  code 
is  documented  and  Intermediate  and  final  results  using  the  developed 
form  factor  are  presented  and  compared  with  other  work  in  the  liter¬ 
ature.  Because  of  continuous  code  developments  at  AFGL.the  non-LTE 
and  COj  continuum  modifications  were  incorporated  into  FASC0D1B* (see 
Appendix  B) ,  and  the  resulting  code  is  designated  in  this  report  FASCOD2-82. 

*The  CO2  lineshape  modifications  developed  in  this  report  are  presently 
available  to  the  public  in  FASC0D1C. 
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Section  5  presents  a  discussion  of  the  layering  procedure 
used  in  FASCOD  and  how  it  may  be  improved.  Such  suggestions  are  not 
incorporated  into  the  code.  Finally,  Section  6  contains  a  few  sum¬ 
mary  remarks  on  how  the  various  modules  and  subroutines  developed  here 
will  be  merged  into  an  integrated  code  called  FASC0D2. 


2.  FASCOD2-80  GEOMETRY  AND  AEROSOL  MODELS 


Here  the  radiative  properties  represented  by  the  line  modeling  in 
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FASCODE  ’  are  merged  with  gaseous  and  aerosol  atmospheric  models  and 

the  new  spherical  geometry  package10  to  produce  a  user  oriented  and 

automated  code  that  calculates  the  atmospheric  transmittance  and  radiance. 

The  user's  guide  to  the  code  and  the  results  from  example  calculations 

are  presented  in  Appendix  A. 

2.1  Merging  ATMPTH  With  FASC0D1A:  A  new  geometry  routine  called  ATMPTH10 
was  developed  at  AFGL  as  an  off-line  program  that  prepared  an  atmospheric 
data  file  for  FASCODE.  The  output  of  ATMPTH  contains  integrated  absorber 
amounts  and  density  weighted  average  temperatures  and  pressures  for  the 
homogeneous  layers  traversed  by  the  optical  path.  Prior  to  the  develop¬ 
ment  of  ATMPTH,  programs  called  LAYER10  and  DRIVER10  were  run  sequentially 
to  prepare  the  same  data.  ATMPTH  is  more  flexible  and  user  oriented  than 
its  predecessors.  A  choice  of  seven  (including  one  that  may  be  supplied 
by  the  user)  atmospheric  models  are  provided,  and  the  user  specifies  the 
placement  of  the  layer  boundaries.  To  automate  execution,  ATMPTH  is  made 
an  on-line  subroutine  of  FASCODE,  and  a  discussion  of  this  process  follows. 

Operation  of  ATMPTH  and  its  subroutines  is  controlled  by  three  input 
cards  whose  format  is  described  in  Appendix  A.  The  first  input  card 
specifies  the  choice  of  model  atmosphere.  There  are  six  built-in  atmos¬ 
pheres  specified  by  pressure,  temperature  and  the  molecular  densities  of 

HjO  CO2,  0^,  ^0,  CO,  CH^,  O2  and  N2  at  34  altitudes.  These  model  at- 
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mospheres  are  essentially  the  same  as  those  in  LOWTRAN  ,  except  that  the 
densities  are  given  in  different  units.  One  of  the  six  model  atmospheres 
stored  in  subroutine  MDLATM,  or  a  non-standard  user-supplied  atmosphere  as 


read  by  subroutine  NSMDL,  is  selected  and  stored  in  the  common  block 
/MDATA/ . 

The  second  control  card  must  contain  three  of  the  following  five 

path  parameters:  HI,  H2,  ANGLE,  RANGE,  and  BETA.  The  reader  should 

9 

refer  to  the  LOWTRAN  report  for  allowable  combinations  of  these  parameters. 
Next,  the  desired  FASCODE  layer  boundaries  are  read  in.  These  need  not 
correspond  to  the  altitudes  contained  in  the  model  atmospheres  or  those 
of  Hi  and  H2,  the  altitudes  of  the  path  end-points.  The  only  restriction 
on  the  boundaries  is  that  they  be  chosen  such  that  the  ratio  of  the  mean 
half-widths  in  adjacent  layers  is  not  greater  than  2  to  1.^  A  new  atmos¬ 
pheric  profile  containing  both  the  previous  altitudes  and  the  desired 
FASCODE  boundaries  is  formed.  Pressures,  temperatures,  and  molecular 
densities  at  these  boundaries  are  found  by  interpolation  of  the  model 
atmosphere. 

Subroutine  GMTRY  is  called  next  to  calculate  the  refracted  path  through 
the  atmosphere  and  to  integrate  the  amounts  of  the  absorbers  along  that  path. 
Subroutine  RFRPTH  drives  the  integration  and  requires  as  input  the  set  Hi, 
ANGLE,  H2  and  LEN.  These  are  determined  by  GMTRY  using  the  three  paths 
parameters  supplied.  A  final  atmospheric  profile,  generated  by  RFRPTH, 
starts  at  HMIN  (tangent  height  if  any)  and  goes  up  to  MAX  (Hi,  H2)  (the 
larger  of  Hi  and  H2).  The  final  profile  is  specified  at  HI,  H2  and  all 
of  the  model  atmosphere  altitudes  and  FASCODE  layer  boundaries  between 
HMIN  and  MAX  (HI,  H2) .  The  determination  of  the  refracted  path,  integration 
of  the  absorber  amounts,  and  formation  of  the  homogeneous  layers  is  done 
using  all  of  the  altitudes  in  the  final  profile  as  layer  boundaries.  Once 
the  path  and  integrated  amounts  are  known,  layers  are  merged  in  ATMPTH  to 
form  the  desired  FASCODE  layers  which  are  the  boundaries  HMIN,  Hi,  H2  and 
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the  FASCODE  boundaries  between  HMIN  and  MAX  (HI ,  H2) .  Note  that  the  calculations 
of  the  refracted  path  and  absorber  amounts  are  done  using  the  model  Atmosphere 
altitudes  as  additional  layer  boundaries.  Only  after  the  calculation  pro¬ 
cedure  is  complete,  is  the  number  of  layers  reduced.  This  is  done  to 
avoid  degrading  the  model  atmosphere  data  in  the  process  of  performing 
the  path  integrals.  A  general  flow  diagram  on  ATMPTH  is  shown  in  Figure  2.1. 

The  addition  of  ATMPTH  as  a  subroutine  of  FASCODE  was  straightforward. 

A  single  statement  in  the  FASCODE  main  routine  was  added  to  call  ATMPTH 
one  time.  During  the  execution  of  this  call  statement,  all  of  the  required 
atmospheric  data  is  generated  and  stored  in  the  common  blocks  /PATH/  and 
/Oi.TPUT/,  which  also  appear  in  FASCODE  subroutine  PATH.  PATH  is  called  as 
before,  first  to  initialize  the  binary  line  absorption  file  and  then  once 
for  each  layer  to  define  the  pressure,  temperature,  absorber  amounts,  and 
half-widths  for  that  layer.  Two  changes  in  PATH  are  worth  noting  at  this 
point.  First,  a  new  section  of  code  was  added  to  create  and  save  the  mean 
half-widths  for  all  of  the  layers.  This  portion  of  code  is  executed  during 
the  first  call  of  PATH  when  the  line  file  is  initialized.  Second,  layer  data 
are  now  defined  in  terms  of  the  stored  array  values.  Statements  of  the  following 
form  have  beer,  added  to  the  section  of  PATH  that  is  executed  for  each  layer 

PAVE  -  Pti.VR  (L)  (average  pressure)  (2.1) 

TAVE  ■  'I BAR  (L)  (average  temperature)  (2.2) 


1. 

2. 


3. 

5. 

6. 

7. 

8. 


9. 


10. 
1  1 . 


INTEGRATED  AMOUNTS  FOR  FASCODl  LAYERS 


FI  (JURE  2.1  FLOW  DIAGRAM  OF  SUBROUTIME  ATMPTH 
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where  L  is  the  layer  index.  These  quantities  were  rreviously  read  from 
a  file  generated  by  an  off-line  run  of  ATMPTH.  The  ATMPTH  subroutine 
has  no  further  interaction  with  FASCODE  beyond  the  assignment  of  the 
average  layer  quantities  in  PATH. 

To  verify  that  ATMPTH  had  been  successfully  Interfaced  with  FASCODE 
as  a  subroutine,  comparison  runs  for  horizontal,  vertical,  and  slant 
(long  and  short)  paths  were  made  using  the  latest  version  of  FASCOD2-80 
with  ATMPTH  as  a  subroutine  and  an  original  version  of  FASC0D1  using  ATMPTH 
as  an  off-line  program.  The  codes  produced  identical  output  in  all  cases, 
indicating  that  the  atmospheric  information  was  being  correctly  transferred 
to  FASCOD2  from  the  ATMPTH  subroutine. 

2.2  Adding  the  LOWTRAN  Aerosol  Models  to  FASC0D1A:  In  its  original  form, 
FASC0D1A  did  not  Include  the  effects  of  molecular  scattering  or  aerosol 
scattering  and  absorption  outside  the  microwave  region.  This  extinction 
mechanism  or  zeroth-order  scattering  has  been  added  as  part  of  the  de¬ 
velopment  of  FASCOD2-80.  A  single  expression  for  the  molecular  scattering 

coefficient  along  with  the  complete  set  of  AFGL  aerosol  models  were  taken 
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directly  from  th"  LOWTRAN  5  computer  code.  Linkage  with  FASC0D1  is  ac¬ 
complished  with  a  newly  written  subroutine,  AERSOL,  which  serves  as  a 
driver  for  LOWTRAN  routines  AERPRF,  PRFDTA,  EXABIN,  EXTDTA,  and  AEREXT. 

The  only  changes  made  to  these  routines  were  to  reduce  the  size  and  number 
of  common  blocks  wherever  LOWTRAN  variables  were  not  needed  by  FASCODE. 
AERSOL,  its  related  subroutines,  and  the  modifications  made  to  FASC0D1A 
to  incorporate  the  new  extinction  mechanism  are  documented  here. 

AERSOL  is  first  called  to  generate  the  aerosol  densities  at  the  34 
standard  altitudes  of  the  model  atmospheres.  During  this  initial  call  a  new 
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aerosol  control  card,  described  In  Appendix  A,  Is  read.  Subroutine  AERPRF 
Is  then  called  34  times  (once  for  each  altitude)  to  load  the  appropriate 
aerosol  density  profile  from  the  data  stored  in  subroutine  PRFDTA  into  the 
array  EHM(34)  in  common  block  /PROF/ .  These  densities  are  unitless  and  are 
used  to  obtain  equivalent  sea  level  absorber  amounts.  At  this  point,  control 
returns  to  .FASCODE 's  main  routine  where  ATMPTH  is  called  to  generate  the  layer 
data.  EHM(34)  is  passed  to  ATMPTH  where  it  is  handled  in  the  same  manner 
aa  the  molecular  density  profiles.  The  integrated  aerosol  amounts  are 
stored  in  the  array  AWKAER(L)  and  are  later  assigned  to  the  simple  variable 
name  WKAER  in  subroutine  PATH.  The  primary  function  of  the  AERSOL  routine, 
namely  to  generate  the  aerosol  density  profile  and  pass  it  to  ATMPTH,  is 
complete  at  this  point. 

AERSOL  is  subsequently  called  once  for  each  layer.  Its  function  now 
is  to  load  the  FASCODE  arrays  ABSRB  and  SCTTR  with  the  additional  required 
optical  depths.  ABSRB  was  previously  used  by  FASCODE  to  store  molecular 
continuum  optical  depths  over  the  wavenumber  range  of  VIABS  to  V2ABS  in 
increments  of  DVABS.  The  aerosol  extinction  and  molecular  scattering 
optical  depths  were  added  to  this  array  because  they,  like  molecular  con¬ 
tinuum  absorption,  are  relatively  slow  functions  of  wavenumber.  The 
user  should  note  that  the  array  name  "ABSRB"  is  now  somewhat  of  a  misnomer 
since  scattering  optical  depths  are  also  included.  Renaming  the  array 
"EXTNCT"  or  "ABSSCT"  to  avoid  any  ambiguity  is  left  to  the  user's  discretion. 
The  array  SCTTR  covering  the  same  wavenumber  region  at  the  same  resolution 
was  empty  prior  to  the  current  addition  of  the  molecular  and  aerosol 
scattering  optical  depths.  This  array  had  been  incorporated  into  an  earlier 
version  of  the  code  with  foresight  of  the  current  modifications. 


Optical  depths  for  the  new  mechanisms  are  defined  as  the  product 
of  the  Integrated  amount  and  the  appropriate  attenuation  coefficient. 

The  aerosol  amount  for  a  given  layer  is  WKAER.  The  coefficients  are 
functions  of  the  aerosol  model,  the  relative  humidity  (for  a  altitudes 
between  0  and  2  Km),  and  the  wavenumber.  During  AERSOL's  first  call  within 
the  FASCODE  layer  loop,  EXABIN  is  called  to  obtain  a  set  of  altitude  and 
wavelength  dependent  attenuation  coefficients  using  the  data  stored  in 
subroutine  EXTDTA.  The  renultlnc  absorption  and  scattering  coefficients 
are  stored  in  arrays  ABSC  and  EXTC.  A  wavenumber  loop  in  AERSOL  calls 
AEREXT  repeatedly  to  determine  the  coefficients  at  the  desired  wavenumbers 
by  interpolation  of  the  ABSC  and  EXTC  arrays.  The  aerosol  extinction  and 
scattering  optical  depths  are  then  defined  by  expressions  of  the  form: 


T.  _  -  EXT  *  WKAER  (2.3) 

A  |  u 

and 

ta  s  -  (EXT  -  ABS)  *  WKAER  .  (2.4) 

If  a  FASCODE  layer  is  not  totally  contained  within  one  of  the  four  aerosol 
altitude  regions  (0-2  Km,  2-10  Km,  10-30  Km,  or  30-100  Km),  a  linear  com¬ 
bination  of  coefficients  is  used.  Molecular  scattering  is  handled  somewhat 
differently,  although  the  basic  idea  is  the  same.  Molecular  amounts  are 
measured  in  terms  of  equivalent  sea  level  distances  defined  as 


WKMOL  -(£)  (f)s. 


(2.5) 
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where  P  and  T  are  the  molecular  density  weighted  average  pressure  and 
temperature  for  the  layer.  The  quantities  P0  and  Tq  refer  to  STP,  and  S 
is  the  geometric  path  length  through  the  layer.  The  scattering  coefficient 
(per  unit  of  sea  level  distance)  is  modelled  by  the  following  empirical 
expression.9 


SCTM  -  v4/ (9. 26799  *  1018  -  1.07123  *  109  *  v2)  . 


(2.6) 


The  molecular  scattering  optical  depth  is  then  given  by  the  product  of  the 
coefficient  and  the  "amount"  as 


T  -  SCTM  *  WKMOL  . 
m,  s 


(2.7) 


Loading  the  ABSRB  and  SCTTR  arrays  for  each  layer  concludes  the  interaction 
of  AERSOL  with  FASCODE'S  main  program.  A  flow  diagram  of  AERSOL  and  its 
subroutines  is  shown  in  Figure  2.2. 

As  a  result  of  adding  the  scattering  mechanisms  to  FASCODE,  the  radiance 
algorithm  required  modification.  The  Rl  array  which  contains  the  optical 
depths  on  the  finest  wavenumber  grid  ultimately  has  the  ABSRB  array  merged 
into  it.  Rl,  therefore,  now  contains  total  extinction  optical  depths 
rather  than  just  absorption.  The  transmission  algorithm  involves  only 
exponentiation  of  the  Rl  array  and  remains  unchanged.  The  radiance  expression 
in  subroutine  EMIN  was  originally  written  as 


BBRAD  (I)  -  (l-TR(I) )  *  BB 


(2.8) 


J 


and  la  now  incorrect  because  TR(1)  Is  an  extinction  rather  than  an 
absorption  transmittance.  The  expression  is  easily  corrected  by  multi¬ 
plying  it  by  the  ratio  of  absorption  to  extinction.  It  turns  out  to 
be  more  convenient  to  use  unity  minus  the  ratio  of  scattering  to 
extinction.  This  is  expressed  in  the  code  as 

BBRAD(I)  -  (l-TR(I))  (1-ALB)  *  BB  ,  (2.9) 

where  ALB  Is  the  ratio  of  scattering  to  extinction  optical  depths,  i.e.  the 
single  scattering  albedo.  Notice  that  equation  (2.9)  reduces  to  equation 
(2.8)  in  the  limit  of  zero  scattering. 

To  verify  that  the  aerosol  routines  were  functioning  correctly  as 

q 

part  of  the  FASCODE  program,  comparisons  were  nu.be  with  LOWTRAN  results. 
FASCODE  and  LOUTRAN  now  contain  identical  aerosol  models,  but  the  molecular 
absorption  processes  are,  of  course,  completely  different.  To  obtain  com¬ 
parable  data  from  FASCODE,  the  molecular  absorption  was  artiflcally  suppressed 
and  additional  output  was  written.  Vertical  paths  were  chosen  to  eliminate 
any  differences  caused  as  a  result  of  the  different  geometry  routines  used 
in  the  two  codes.  Comparisons  of  aerosol  absorber  amounts,  coefficients, 
and  transmlttances  indicated  excellent  agreement  between  FASCODE  and  LOWTRAN. 

A  slight  different  in  the  aerosol  coefficient  in  the  0-2  Km  boundary  layer 
region  wsb  observed  and  attributed  to  the  different  relative  humidity  cal¬ 
culations  within  the  two  codes.  The  relative  humidity  algorithm  developed 
under  this  contract  for  use  with  FASCODE  differs  from  the  one  currently  used  in 
LOWTRAN.  Although  the  net  effect  on  the  aerosol  coefficients  is  small ,  typically 
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about  one  percent,  it  does  lead  to  differences  in  the  attenuation  coefficients 
when  used  as  a  reference  point  for  interpolation  of  the  stored  data.  The 
following  brief  description  of  the  two  algorithms  will  illustrate  the 
fundamental  differences  and  how  they  may  be  eliminated. 

The  LOWTRAN  geometry  routines  integrate  the  relative  humidity,  weighted 
by  the  aerosol  density,  while  performing  the  usual  path  integrations.  The 
calculation  is  the  numerical  approximation  of 


RH 


RH(h)  *  WK(h)  dZ 


WK(h)  dZ 


(2.10) 


where  RH  and  WK  are  arrays  containing  the  relative  humidity  and  aerosol  density 
profiles  as  a  function  of  altitude  h,  and  dZ  is  the  incremental  path  length.  Both  the 
numerator  and  denominator  ofEq.  (2.10)  are  formed  during  the  geometry  cal¬ 
culations  and  then  their  quotient  is  formed.  By  contrast  the  FASCODE  relative 
humidity  calcdlat ion  occurs  after  all  geometry  calculations  are  completed.  A  discrete 
sum  of  layer  properties  rather  than  "integration"  of  atmospheric  profiles 
is  uBed.  The  average  relative  humidity  is  expressed  as 


RH(L)  WK(L) 

-  ,  (2.11) 

2j  WK(L) 

0-2  Km 


RH 


0-2  Km 
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where  RH  and  WK  are  arrays  containing  the  relative  humidity  and  aerosol 
absorber  amounts  for  the  FASCODE  layers.  Equations  (2.10)  and  (2.11) 
both  represent  an  aerosol  weighted  relative  humidity.  As  it  turns  out, 
the  denominators  of  the  two  expressions  are  numerically  equivalent.  The 
difference  lies  in  the  numerators,  with  the  sum  in  Equation  (2.11)  being 
an  approximation  of  the  integral  in  Equation  (2.10).  Although  time  and 
resources  did  not  permit  the  replacement  of  FASCODE 'S  relative  humidity 
algorithm  with  the  more  accurate  one  found  in  LOWTRAN,  a  method  to  do  so  is 
outlined  below  for  future  reference: 

1.  Provide  additional  array  space  for  the  relative  humidity  profile  and 
a  weighted  relative  humidity  sum. 

2.  Add  a  section  of  code  prior  to  the  geometry  calculations  to  generate 
and  load  the  relative  humidity  profile. 

3.  Modify  the  geometry  routines  to  perform  an  additional  integration, 
namely,  the  relative  humidity  weighted  by  the  aerosol  density. 

Note,  this  integration  need  only  cover  the  0-2  Km  boundary  layer 
region. 

4.  After  the  geometry  calculations  are  completed,  divide  the  weighted 
relative  humidity  sum  by  the  boundary  layer  aerosol  amount. 


3.  NON-LTE  RADIANCE  AND  TRANSMITTANCE:  FASCOD2-82 


In  this  section  we  review  the  monochromatic  radiative  transfer  proper¬ 
ties  of  an  absorbing  gaseous  medium  which  may  not  be  in  local  thermodynamic 
equilibrium  (LTE).  Those  assumptions  and  relationships  which  are  valid 
only  under  LTE  circumstances  are  distinguished  from  the  quite  general 
aspects  of  radiative  transfer.  The  transfer  equation  describes  the  influ¬ 
ence  which  the  atmospheric  environment  (the  optical  medium)  has  on  the 
radiation  field  through  both  losses  and  emission,  the  problem  which  FASC0D2 
has  been  designed  to  address.  The  radiation  field  also  affects,  in  various 
ways,  the  state  of  the  gas,  both  as  a  driving  force  that  stimulates  molecu¬ 
lar  transitions  and  as  a  mechanism  for  energy  exchange  and  relaxation 
within  the  medium.  The  largest  single  effect  of  solar  radiation  on  the 
earth's  atmosphere  is  simple  thermal  heating.  However,  in  the  upper  atmos¬ 
phere  selectively  enhancM  molecular  excitations  and  sun-driven  photochemi¬ 
cal  processes  combine  with  the  slow*:  relaxation  rates  of  the  rarefied  gas 
medium  to  drive  the  medium  away  from  3  state  of  thermodynamic  equilibrium. 

3.1.  Theoretical  Development:  High  altitude  radiance  measure¬ 

ments  have  confirmed  that  certain  species,  notably  was.  ?r  "apor,  carbon 
dioxide,  ozone,  end  nitric  oxide,  may  be  found  in  non-LTE  states  above  an 
altitude  of  about  40  -  60  km  (even  lower  for  NO).^  In  the  upper 
atmosphere,  it  is  necessary  to  account  for  these  non-equilibrium  states  of 
the  optical  medium  if  the  radiative  transfer  is  to  be  properly  modeled.  At 
lower  altitudes,  the  I.TE  assumption  is  generally  valid,  although  the  local 


temperature  remains  a  strong  function  of  altitude  from  the  boundary  mixing 
layer  to  the  upper  stratosphere. 

FASC0D2  calculates  only  the  unscattered  radiative  intensity  for  a 
particular  direction  of  observation  using  known  atmospheric  properties. 

The  code  therefore  implicitly  decouples  the  problem  of  solving  for  these 
intensities  from  the  problem  of  finding  the  state  of  the  medium  itself. 

For  conditions  of  local  thermodynamic  equilibrium,  the  properties  of  an 
atmosphere  which  are  needed  to  perform  these  calculations  (using  LTE  rela¬ 
tionships)  are  accessable  in  the  form  of  stored  atmospheric  profiles  (tem¬ 
perature,  pressure,  gas  densities)  plus  an  off-line  file  containing  mole¬ 
cular  line  data  at  T“2968K*.  The  non-LTE  problem  however,  requires  a 
more  detailed  description  of  the  propagating  medium,  and  removal  of  all  LTE 
assumptions  used  in  the  radiative  transfer  algorithm. 

The  non-LTE  transfer  problem  can  be  described  by  the  following  quite 
general  governing  equation  for  the  monochromatic  radiative  intensity  Iv: 


K. 


I  + 

V 


[  3.1  ] 


The  variable  s  is  taken  to  be  the  length  along  a  (possibly  curved)  line- 
of -sight  through  the  atmosphere.  (This  could  equally  well  have  been 
defined  a  molecular  amo  in4,  per  unit  cross-sectional  area  or  gas  column 
dens  r>;  )  The  "ocal  monochromatic  volumetric  extinction  coefficient  is 
Kv  ,  ani  represents  the  source  U  .  for  radiative  emission  per  unit 

length.  It  is  important  to  note  that  both  ana  R,,  are  linear  func- 
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tions  of  the  densities  of  each  IR-active  gas  and  aerosol  constituent. 
Therefore,  the  separate  contributions  to  both  quantities  can  be  added 
line-by-line,  molecule-by-molecule ,  aerosol-by-aerosol,  etc. 

The  basic  objective  is  to  calculate  the  contribution  to  both  Kv  and 
Rv  for  all  molecular  absorption/emission  lines.  In  order  to  illustrate 
the  processes  involved,  a  representative  vibrational-rotational  transition 
is  considered  between  an  upper  state  u  and  lower  state  &  corresponding 
to  a  spectral  line  centered  at  frequency  VQ  .  In  terms  of  the  standard 
Einstein  coefficients  B^u  ,  Bu^  ,  and  ,  for  absorption,  stimu¬ 

lated,  and  spontaneous  emission,  respectively,  the  volumetric  absorption 
and  radiance  for  this  transition  are  given  by 


■  S  «*,  V  *v  -  V  1  3-2 

and 

hv  13.3 

R  ■  A„  n  X  • 

V  ATT  Hu  U  V 

The  quantities  ^  ,  4^  ,  and  Xv  represent  distinct  lineshepe  functions 

of  the  spectral  distance  (  v  -  v  ) ,  and  are  referenced  to  one  another 
by  having  4>  (0)  =  (0)  =  X  (°)  •  Number  density  populations  of  the 

lower  and  upper  V-R  states  are  given  by  n.  and  n  ,  respectively. 

Ay  U 

The  Einstein  coefficients  are  a  set  of  rate  constants  which  are  sped 
fic  to  a  particular  transition,  but  which  are  not  independent.  This  fact 
is  generally  true  and  most  easily  demonstrated  by  considering  a  situation 
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of  complete  equilibrium.  At  equilibrium  the  radiation  field  within  a 
sufficiently  large  medium  of  temperature  T  is  isotropic  and  has  a  spec 
tral  intensity  given  by  the  Planck  function  B^CT)  in  the  form 


Bv  (T) 


2 

c 


,  hv/kT 
(e 


-1)  • 


[  3.4  ] 


Additionally,  the  populations  of  the  upper  and  lower  states  have  the  ratio 


nu 


-hv  /k'T 

O 


l  3.5  ] 


given  by  the  Boltzmann  distribution  for  the  temperature  T.  The  constant 
integer  g's  are  single  3tate  degeneracy  factors.  The  superscript  "e"  will 


be  used  throughout  to  denote  the  equilibrium  state  for  a  specified  tempera¬ 
ture.  Since  the  intensity  doesn't  change  with  position  in  this  case,  the 


equation  of  transfer  reduces  to 


\  •  K  \  (T>  • 


I  3.6  ] 


which  implies  that  at  equilibrium 


[  3.7  ] 


Furthermore,  the  fundamental  symmetry  of  the.  quantum  matrix  elements  for 

the  processes  of  absorption  and  stimulated  emission  requires  that  B  = 

g  B  .  Considering  Equations  3.4,  3.5,  and  3.7  evaluated  at  the  line 
l  Hu 

center  point  V  *  VQ,  the  three  Einstein  coefficients  can  be  expressed  as 
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where  B  is  a  single  rate  coefficient. 


The  above  argument  can  also  be  used  to  relate  the  lineshapes  to  one 

another.  There  is  a  primary  relationship  between  emission  lineshapes 

Xv  =  which  follows  from  the  Wiener-Khintchine  theorem  and  possibly 

13  14 

more  general  arguments  as  well.  1  One  also  has  from  Equations  3.5» 
3.7,  and  3 . 8  that 


therefore  requiring  only  a  single  lineshape  model  to  be  used  in  describing 
the  three  radiative  processes.  Arguments  similar  to  the  one  above  have 
been  made  previously.  Gilles^  develops  the  corresponding  lineshape 
relationships  using  a  different  lineshape  normalization  condition  and 
somewhat  different  notation  throughout. 


In  summary,  the  equilibrium  single-line  radiative  problem  within  a 
given  medium  can  be  adequately  described  in  terms  of  a  single  rate  coef¬ 
ficient  (or  linestrength)  and  a  single  lineshape  assumed  known  from  other 
considerations.  The  above  Einstein  coefficients  can  be  applied  to  the 
non-LTE  problem  with  equal  validity  since  they  represent  molecular  proper¬ 
ties.  Any  one  of  several  assumptions  can  be  made  regarding  how  the  line- 
shapes  should  be  related  to  each  other  under  non-LTE  conditions,  that  is, 
how  (3.10)  can  be  generalized.  In  view  of  the  fact  that  both  collisional 
and  Doppler  line  broadening  mechanisms  depend  largely  on  the  bulk  thermo¬ 
dynamic  properties  of  pressure  and  temperature,  and  not  on  other  degrees  of 
freedom  such  as  the  particular  V-R  state  populations  n.  and  n  ,  it 
seems  reasonable  to  assume  that  the  lineshapes  should  have  their  LTE  form 
even  under  non-LTE  conditions.  One  then  only  requires  detailed  knowledge 
of  all  V-R  single  state  populations  to  solve  the  radiative  tranfer  problem 
by  applying  Equations  3.2,  3.3,  3.8,  and  3.10. 

In  carrying  out  calculations  based  on  the  above  theory,  it  becomes 
rather  awkward  to  actually  allow  and  to  be  different  lineshape 

functions.  The  overall  (true  absorption  iiinus  induced  emission)  lineshape 
for  absorption  is  a  weighted  difference  between  the  functions  4>v  and 
\py  ,  and  would  therefore  change  as  the  upper  to  lower  state  population 
ratio  changes  with  altitude  and  from  line  to  line.  To  avoid  this  compli¬ 
cation,  the  approximation  is  made  that  the  overall  (net)  absorption  line- 
shape  is  the  same  under  rton-LTE  conditions  as  it  is  under  LTE 
conditions.  From  (3.10)  it  can  be  seen  that  this  approximation  introduces 
errors  in  the  lineshape  wings  of  order  exp[-h(  v  -  VQ  )/kT).  Except  for 
the  microwave  region,  the  errors  introduced  by  this  lineshape  approximation 


are  quite  small  in  comparison  to  the  typical  uncertainties  in  the  non-LTE 
1  inestrengths ,  The  differences  are  especially  insignificant  for  the  narrow 
Doppler  lines  found  at  high  altitudes. 


With  the  above  approximation,  the  single-line  volumetric  absorption 
and  local  thermal  radiance  is  given  by 


hv  -  ,  v  . 

r B  <BuV  H  nu>  *v 


hV  s 

^2  B 


[  3.11  1 


[  3.12  ] 


This  result  is  a  general  one  for  a  two- level  system.  The  lineshapes  are 
functions  of  the  bulk  thermodynamic  quantities  such  as  pressure, 
temperature,  plus  individual  gas  densities  which  can  separately  influence 
the  far  wing  lineshapes.  The  populations  are  given  by  the  Boltzman 
distribution  only  under  conditions  of  local  thermodynamic  equilibrium.  In 
the  presense  of  a  strong  radiation  field,  or  under  conditions  where 
relaxation  rates  are  small,  the  upper  and  lower  state  populations  will 
often  bear  no  simple  relationship  to  each  other.  Therefore,  the  absorption 
and  radiance  owing  to  any  single  transition  (and  to  all  transitions  of  all 
molecules)  are  no  longer  simply  related  for  a  non-LTE  atmosphere,  and  more 
specifically,  the  radiance  from  such  an  optically  thick  atmosphere  will  not 
be  given  by  the  Planck  function. 


In  the  upper  atmosphere,  relaxation  rates  for  the  different  excita¬ 
tion  modes  of  various  molecules  are  found  to  vary  considerably.  Rotational 


degrees  of  freedom  generally  relax  much  more  rapidly  to  a  state  of  thermal 
equilibrium  than  do  vibrational  modes.  The  work  of  Degges  and  Smith^ 
makes  use  of  this  fact  in  order  to  construct  a  high  altitude  model 
which  calculates  vibrational  state  populations,  while  assuming  that  the 
rotational  state  populations  are  thermally  distributed  within  each  vibra¬ 
tional  state  according  to  the  Boltzman  distribution  for  the  translational 
(bulk)  temperature.  This  model  employs  a  single  kinetic  temperature  which 
is  a  function  of  altitude,  plus  a  set  of  vibrational  state  population 
variables  which  are  found  by  solving  a  model  for  the  coupled  gas-radiation 
field  problem. 

These  vibrational 'population  variables  can  be  used  to  relate  the 
non-LTE  single-line  absorption  and  radiance  to  the  same  quantities  calcu¬ 
lated  for  an  atmosphere  in  equilibrium  at  the  kinetic  temperature.  The 
relationship  between  equilibrium  and  non-equilibrium  radiative  properties 
for  a  single  line  can  be  summarized  by 


and 


8u  nfc  ~  8£  nu 

e  e 

8u  n£  -  8a  nu 


V 


[  3.13  ] 


R  -  R  e 

Rv  e  Rv 

n 

u 
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which  can  also  be  written  in  terms  of  the  Planck  function  as 


Rv  -  \  vT)  V 

n 

u 


[  3.15 
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In  the  upper  atmosphere  under  non-LTE  conditions,  some  vibrational 
state  populations  are  increased  by  as  much  as  a  factor  of  100  to  1000  over 
LTE  populations.  The  ground  (lowest  energy)  vibrational  state  is  usually 
only  marginally  depleted.  The  net  effect,  as  can  be  seen  from  3.13  and 
3.15,  tends  to  be  a  slightly  diminished  spectral  absorption  coefficient, 
while  the  radiance  tends  to  be  markedly,  though  selectively,  enhanced  in 
some  bands  of  ozone,  carbon  dioxide,  nitric  oxide,  and  water  vapor.  The 
other  major  gases  are  also  affected  to  a  lesser  extent. 

Some  sample  calculations  have  been  made  using  this  high  resolution 
non-LTE  model.  Figure  3.1  shows  the  enhancement  of  the  predicted  radiance 
by  non -equilibrium  vibrational  populations  which  would  be  seen  by  an 
observer  looking  along  a  500  km  horizontal  path  at  an  altitude  of  100  km, 
based  on  the  1962  U.S.  Standard  Atmosphere  in  the  model  of  Degges  and  Smith.  The 
comparison  of  LTE  and  non-LTE  calculations  is  made  for  a  spectral  region 
within  the  4.3  micron  band  cf  CO^.  This  band  and  others  often  show  radi¬ 
ance  enhancements  by  several  orders  of  magnitude  over  what  could  be  expec¬ 
ted  of  a  gas  at  thermal  equilibrium. 


3.2.  Non-LTE  Code  Modifications:  The  relations  derived  in  the  previous 
section  provide  for  the  calculation  of  absorption  and  radiance  spectral 
functions  using  two  separate  line-by-line  convolutions.  (The  two  sums 
would  be  related  by  the  Planck  function  in  the  equilibrium  limit.)  The 
non-LTE  algorithm  of  FASC0D2  calculates  Kv  directly,  but  the  radiance  is 
computed  from  a  line-by-line  sum  of  the  quantity  Cy  defined  by 
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Figure  3*1  Radiance  in  Units  of  Eauivilent  Blackbody 
Temperature  for  a  500  km  Path  at  100  km  Altitude  Assuming 
LTE  and  Non-LTE  Vibrational  Populations 
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[  3.16  ] 


r  rn  V  -  ■■ 

Cv  Kv  BV(T) 

i. 

which  can  be  viewed  as  representing  that  part  of  the  non-LTE  absorption 
which  is  not  returned  to  the  radiation  field  as  thermal  emission.  Any  two 
of  the  spectral  functions  ,  R^  ,  and  can  easily  be  used  to 

obtain  the  third,  so  the  choice  of  summing  rather  than  has  been 

made  for  computational  reasons.  can  be  seen  to  vanish  in  the  equili¬ 

brium  limit,  and  furthermore,  in  an  atmosphere  of  mixed  gases,  only  transi¬ 
tions  of  those  gases  which  are  not  in  thermal  equilibrium  will  contribute 
to  This  property  is  used  by  FASC0D2  to  shorten  the  second  line-by- 

line  sum  where  possible, 


To  compute  the  absorption  and  radiance  for  optical  path  segments  in 
the  upper  atmosphere,  FASC0D2  employs  the  vibrational  population  model  of 
Degges  and  Smith  ^ .  Each  IR  line  corresponds  to  a  transition  be¬ 
tween  two  molecular  vibrational  states.  The  identification  of  these  two 
states  is  made  by  program  BCDNLTE,  a  separate  routine  modeled  after  BCDMRG 
which  prepares  a  blocked  line  data  file  for  FASC0D2.  The  popula¬ 
tion  of  each  vibrational  state  is  referenced  to  its  equilibrium  value, 
calculated  at  the  kinetic  temperature  T.  The  ratios  of  non-LTE  vibrational 
populations  to  their  reference  values  are  given  by  the  enhancement  (or 
depletion)  ratios 


r  _ 


and 


[  3.17  ] 


The  fact  that  rotational  sub-states  are  thermally  distributed  has  been  used 
to  relate  the  vibrational  state  population  ratios  directly  t.o  single  V-R 
state  population  ratios.  That  is,  rotational  partition  functions  and 
Boltzman  factors  produce  a  factor  of  unity  in  the  V-R  ratios.  The  single- 
line  contributions  to  and  can  be  expressed  simply  in  terms  of 

these  ratios  according  to 


and 


where  A  is  the  Boltzman  factor, 
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[  3.18  ] 


t  3.19  ] 


[  3.20  ] 


The  two  factors  in  brackets  can  be  viewed  as  non-LTE  lineatrength  correc¬ 
tion  factors  for  these  absorption  and  "effective  absorption"  line-by-line 
calculations . 


To  summarize,  if  "i"  denotes  a  single  absorption  line  of  a  particular 
molecular  species,  the  two  functions  Kv  and  are  obtained  by  per¬ 

forming  the  two  sums 

Ku  and  cv  ‘2  C  1  3  21 
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which  are  taken  over  all  lines  of  all  molecules  in  the  case  of  ,  and 

over  all  lines  connecting  non-equil ibrium  vibrational  populations  in  the 
case  of  Cv  .  Both  Kv  and  C^,  employ  the  same  approximate  equilibrium 
lineshape  functions  as  discussed  above.  The  complete  solution  of  the 
radiative  transfer  equation  for  a  path  of  length  L  through  each  homo¬ 
geneous  layer,  including  self -absorption  of  emitted  radiation,  is  given  by 

'v  -  C*'V  +  Bv[l-^][1-*'V]  •  [  3-“  1 

The  effect  of  solving  the  non-equilibrium  problem  in  each  homogeneous  layer 
has  been  to  modify  the  linestrengths  contributing  to  the  local  absorption 
coefficient  ,  and  to  replace  the  Planck  function  Bv  by  the  high 
resolution  function  [  1  -  Cy  /  J . 

In  the  remainder  of  this  section,  the  new  procedures  that  are  neces¬ 
sary  to  make  FASC0D2  perform  the  non-LTE  calculations  are  outlined.  Most 
of  the  modifications  when  considered  individually  have  been  relatively 
minor,  with  the  exceptions  that  one  now  subroutine  was  introduced  (RDPOPS) 
and  an  entire  new  branch  of  the  program  (HIRACQ  and  subroutines)  was  cre¬ 
ated  to  mimic  HIRAC  in  performing  two  line-by-line  calculations  rather  than 
one.  These  modifications  are  described  in  detail  in  Appendix  B.  A  list  of 
the  new  procedures  follows: 

(1)  Identify  and  label  those  upper  and  lower  state  vibrational  states  of 
each  spectral  line  which  are  non-tharmal ly  populated  [BCDNLTE). 
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(2)  Read  the  non-LTE  vibrational  state  population  data  given  by  the  model 
of  Degges  and  Smith,  and  compute  the  population  enhancement  ratio  for  each 
vibrational  state  [RDPOPS], 

(3)  In  computing  Kv  ,  modify  the  single-line  absorption  linestrength 
according  to  Equation  3.18  [HIRACQj. 

(4)  Perform  a  second  line-by-line  calculation  of  the  "effective  absorp¬ 
tion"  function  Cv  patterned  after  the  actual  absorption,  but  based  on 
(3.19),  including  only  those  lines  identified  in  (1)  as  having  either  upper 
or  lower  state  non-equilibirum  populations  [HIRACQ], 

(5)  Modify  the  radiance  calculation  in  each  layer  to  use  the  non-LTE  form 
of  the  Planck  function  Bv  [  1  -  Cv  /  Kv  ][EMIN]. 
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4.  C02  CONTINUUM  ABSORPTION  AND  LINESHAPE  MODELS 

4 

The  line-by-line  HIRACC  molecular  absorption  algorithm  has 
been  extended  to  include  the  lineshape  region  within  256  halfwidths  of  the 
line  center  of  each  spectral  line.^*  ^  The  primary  consideration  in 
extending  the  lineshape  bound  outward  from  64  halfwidths  was  to  better  fit 
the  band  edge  and  region  of  continuum  absorption  thought  to  arise  primarily 
from  the  cumulative  contribution  from  the  far  wings  of  many  overlapping 
lines.  The  continuum  absorption  is  known  to  scale  nearly  linearly  with 
both  the  gas  amount  and  the  broadening  gas  density  (more  precisely,  the 
molecular  collision  frequency)  in  the  high  pressure  (Lorentz  shape)  limit*^ 
It  can  therefore  be  modelled  by  scaling  a  precomputed  or 
measured  continuum  function  based  on  a  reference  pressure  and  temperature. 
The  process  can  be  complicated  if  multiple  broadening  mechanisms  must  be 
accounted  for,  as  in  the  case  of  self  and  foreign  gas  broadening  of  water 
vapor  lines.  The  water  vapor  continuum  has  been  the  focus  of  considerable 
work  at  AFGL  and  is  discussed  in  Reference  [6]. 

In  order  to  model  the  CO^  continuum  at  typical  atmospheric  mixing 
ratios,  it  has  been  assumed  that  CO^  self  broadening  can  be  neglected  in 
favor  of  the  much  stronger  foreign  gas  (principally  nitrogen)  broadening 
effects.  A  number  of  measurements  have  established  that  the  collisional 
broadening  mechanism  (whether  self  or  foreign)  gives  rise  to  a  C02  line 
shape  which  falls  off  much  more  rapidly  in  the  far  wings  than  can  be 
accounted  for  by  the  Lorentz  lineshape  function.^’  ^  This  sub-Lorentzian 


behavior  is  particularly  evident  just  above  the  band  head  of 

the  CO^  4.3  micron  band.  Estimates  of  the  spectral  absorption  in  this 

region,  if  based  on  a  Lorentz  lineshape,  can  be  in  error  by  as  much  as  a 

factor  of  3.  More  recent  atmospheric  transmittance  measurements  made  at 

17  18 

sea  level  and  also  from  high  altitude  balloons  confirm  the 

sub-Lorentslan  behavior  of  line  wings  above  the  4.3  micron  band  edge 

from  about  2380  to  2420  CM-1. 

Several  explanations  for  the  sub-Lorentzian  shape  in  the  line  wings 

have  been  suggested.  They  range  from  arguments  for  a  modified  impact 
19 

approximation  which  is  needed  to  account  for  the  finite  inter¬ 
action  times  of  molecular  collisions,  to  the  suggestion  that  COg -Nj  dimer 

20 

effects  play  an  important  role.  Other  theoretical  models  have 
been  suggested  but  will  not  be  summarized  here.  See,  for  example,  [22] 
and  [23].  These  models  provide  an  important  framework  within  which  the 
measurements  obtained  to  date  can  be  considered,  but  compelling  evidence 
for  any  single  model  has  been  lacking  and  each  model  has  had  limited 
usefulness  as  a  predictive  tool. 

The  most  useful  sub-Lorentzian  lineshape  models  to  date  have  been 
empirical.  The  usual  approach  has  been  to  modify  the  Lorentzian  lineshape 
function  with  a  multiparameter  form  factor  which  is  obtained  by  fitting 
experimental  data  with  trial  synthetic  absorption  spectra.  Winters, 
Silverman,  and  Benedict  (hereafter  shortened  to  WSB)  first  used  an  expo¬ 
nential  form  factor  for  nitrogen  broadening  given  by 


exp 


-  0.46  (| 


v 


5.0) 


0.46 


1 


[  4.1  1 
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which  falls  off  smoothly  starting  at  5  wavenumbers  from  the  line  center 
and  has  two  adjustable  parameters . ^  Susskind  and  Mo  have  also  used 
this  functional  form,  but  they  adjusted  it  to  fall  off  at  only  0.5  wave- 
numbers  from  line  center,  and  fitted  the  exponential  WSB  form  with  somewhat 
different  parameters  using  the  form  factor 

T  0.434-1 

exp  [-  0.601  (|v-v0|  -  0.5)  I  ,  4-J  , 

23 

to  accomplish  the  same  task.  Burch  et.  al.  suggested  that  a  more 

complicated  numerical  form  factor  which  also  starts  to  fall  off  within  one 

wavenumber  of  the  line  center  was  needed  and  they  proceeded  to  construct 

separate  form  factors  for  both  self  and  foreign  gas  broadening  near  the 

three  CO2  bands  at  4.3,  2.6,  and  1.4  microns. ^  Recently,  Cann 

et  al  have  restudied  the  problem  and  suggested  a  different  numerical  form 

24  25 

factor  based  on  similar  considerations.  ’ 

The  available  laboratory  data  are  not  clearly  inconsistent  with  any  of 
these  models,  yet  there  are  seemingly  significant  differences  among  the 
four  form  fa. tors  mentioned  above.  In  particular,  the  shape  of  the  form 
factor  for  >;he  spectral  region  within  10  wavenumbers  of  line  center  varies 
considerably  among  these  models.  Experience  in  building  synthetic  contin¬ 
uum  Absorption  models  such  as  these  has  shown  that  single-line  features  are 
usu*ly  blurred  when  the  wing  contributions  from  many  lines  are  superposed, 
r'iis  fact  suggests  that  it  may  not  be  possible  to  extract  a  unique  form 
factor  from  the  relatively  few  laboratory  band-edge  absorption  measurements 


which  are  available. 


Atmospheric  transmittance  measurements  also  shed  some  light  on  the 
distinctions  among  the  model  predictions.  Smith  et  al  have  recently 
compared  the  continuum  transmission  predicted  by  these  models  to  balloon 
measurements  of  stratospheric  solar  transmission  .  Their  compar¬ 
ison  favors  the  two  models  of  Burch  et  al  and  Susskind  and  Mo,  but  some 

uncertainty  is  introduced  by  the  presense  of  the  pressure  Induced  nitrogen 
absorption  continuum  and  by  the  fact  that  the  strongest  part  of  the  C02 
continuum  above  the  2397  CM-1  band  head  is  compressed  into  a  smaller  spec- 
tral  region  at  the  lower  stratospheric  pressures  where  these  data  are 
taken.  NRL  atmospheric  measurements  taken  over  a  5.1  km  path  at  sea  level 
include  the  same  nitrogen  continuum,  plus  a  somewhat  uncertain  water  vapor 
continuum  contribution  in  the  same  region  and  significant  aerosol  extinc- 
tion  along  the  optical  path.17,  24‘  25  However,  in  spite  of  the 
numerous  overlapping  contributions,  Cann  ot  al  do  claim  an  improved  fit 
to  the  NRL  data  using  their  own  C02  form  factor. 

In  view  of  the  variety  of  empirical  lineshapes  that  have  been  used  to 

fit  the  C02  continuum,  it  was  decided  that  a  new  two-parameter  exponential- 
type  form  factor  particularly  suited  to  the  HIRACC  algorithm  would  be 
studied  for  incorporation  into  FASC0D2.  We  will  briefly  review  below  the 
relevant  features  of  the  HIRACC  line-by-line  calculation  before  detailing 
the  C02  form  factor  which  has  been  used.  The  reader  is  referred  to  [5] 
and  [6]  for  more  complete  descriptions  of  the  HIRACC  algorithm. 

HIRACC  approximates  the  Voigt  lineshape  by  a  weighted  sum  of  Lorentz 
and  Gaussian  (Doppler)  lineshapes.  The  relative  weight  of  each  component 
is  determined  b>  the  Voigt  parameter  <5  ,  defined  in  terms  of  the  Lorentz 


33 


and  Doppler  halfwidths  according  to 


,  “L 


[  4.3  ] 


The  Lorentz  function  is  numerically  constructed  by  the  superposition  of  as 
many  as  four  subfunctions.  A  summary  of  the  functional  forms  for  the 
different  domains  of  each  subfunction  is  given  in  Table  4.1.  Note  that  the 
parameter  X  is  the  distance  from  line  center  measured  in  halfwidths. 

The  notable  functional  forms  are: 

L(X)  Lorentz  function. 


L(X) 

Q4(X) 

Q16(X) 

Q64(X) 

Q256(X) 


The  quartic  function  whose  value,  slope,  and 
curvature  match  the  Lorentzian  at  X"4. 

The  quartic  function  fitting  L(X)  at  X*16. 

The  quartic  function  fitting  L(X)  at  X*64. 

The  quartic  function  fitting  L(X)  at  X«256. 


TABLE  4.1  FASC0D2  LORENTZ  LINESHAPE  SUBFUNCTIONS 
DOMAIN  FI  F2  F3  F4 


4<X<16 

16<X<&4 

64<X<256 


REMAINDER 


qi6-Q64 

Q64-L(X»256) 

L(X»256) 

Q16-Q64 

Q64-L(X«256) 

L(X=256) 

L-Q64 

Q64-L(X“256) 

L(X“256) 

0 

L-L(X*=256) 

L(X=256) 

0 

0 

L 
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The  sum  of  the  four  subfunctions  is  the  Lorentz  function  minus  its  value  at 


X=256  for  the  domain  within  256  halfwidths  of  the  line  center,  while  the 
sum  is  zero  outside  256  halfwidths.  The  column  labelled  "REMAINDER"  shows 
the  "fifth"  function  which  would  be  needed  to  complete  the  Lorentz  shape  in 
each  domain. 

The  advantages  of  the  multifunction  construction  have  been  discussed 
previously"*,  In  particular,  as  the  domain  of  each  successive  sub¬ 
function  grows  larger,  the  functions  FI,  F2,  etc.  grow  successively 
smoother,  so  that  the  wavenumber  grid  on  which  each  function  is  specified 
can  contain  approximately  the  same  number  of  points.  FASC0D2  uses  these 
subfunctions  to  synthesize  the  single  layer  absorptive  optical  depth 
through  line-by-line  summations.  Note  that  each  subfunction  has  a  differ¬ 
ent  characteristic  resolution,  so  that  it  is  advantageous  to  carry  out  five 
separate  line-by-line  sums  in  order  to  construct  five  partial-absorption 
spectra  (designated  Rl,  R2 ,  R3,  R4,  and  ABS  in  the  program)  which  are 
defined  on  successively  coarser  wavenumber  grids.  Once  the  five  convolu¬ 
tions  are  completed,  they  are  merged  into  a  single  optical  depth  function 
which  has  the  high  resolution  characteristic  of  FI. 

The  manner  in  which  the  five  sums  are  actually  computed  deserves  some 
discussion.  The  sums  Rl,  R2 ,  and  R3,  which  make  up  nearly  all  of  the  line 
contribution  to  band  absorption,  are  computed  line-by-line  for  each  absor¬ 
bing  gas  by  HIRACC  and  Its  subroutines.  The  rather  low  resolution  function 
R4  Is  computed  by  subroutine  LBLF4  and  its  subroutines.  These  four  line- 
by-line  convolutions  use  individual  line  strengths  and  halfwidths  which 
are  computed  in  each  atmospheric  layer  from  the  average  layer  pressure  and 
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temperature,  and  the  line  file  data.  The  REMAINDER  or  "fifth" 
subfunction,  however,  has  been  used  only  to  precompute  the  total  contri¬ 
bution  of  all  lines  to  the  continuum  at  standard  temperature  and  pressure. 
When  the  fifth  partial-absorption  function  "ABS"  is  required  in  each  layer, 
this  stored  continuum  function  is  scaled  according  to 

0.5 

ABS(P,T,GNU)  *  ABS(P0,T0,GNU)  [P/PO]  [TO/T]  ,  [  4.4  ] 

which  approximately  reflects  the  dependence  of  the  far  wing  absorption  on 
collision  freque.  „y , 

The  problem  of  finding  an  adequate  model  for  the  far  wing  CO2  line- 
shape  was  first  studied  using  the  simple  exponential  form  factor  (following 
references  [16]  and  [24])  with  the  lineshape 

,  B 

L(X)  EXP[  -  A  (DELTA-C)  ]  for  DELTA  >  C 

L(X)  DELTA  <  C 

having  adjustable  parameters  A,  B,  and  variable  "turn-on"  point  C.  DELTA  is 
the  spectral  distance  from  the  line  center  measured  in  wavenumbers.  Trial 
CO^  absorption  spectra  built  using  the  exponential-type  lineshape  were 
found  to  fit  the  data  progressively  better  as  C  was  decreased  from  the  VSB 
value  of  5  wavenumbers  to  1  wavenumber,  and  finally  to  0.5  wavenumbers. 

The  fits  suggest  that  a  very  small  value  of  C  works  best.  jwever,  one 
knows  that  the  continuum  absorption  is  quite  insensitive  to  the  lineshape 
within  a  few  halfwidths  of  the  line  center.  One  might  conclude  that  the 
"turn-on  point"  should  not  be  chosen  simply  from  considerations  based  on 
fitting  the  continuum,  since  the  small  value  of  C  may  prove  to  be  an  arti- 
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fact  of  forcing  the  form  factor  to  be  of  the  exponential  type  rather  than  a 
meaningful  result.  Additionally,  one  risks  altering  the  integrated  line- 
strength  if  the  "turn-on  point"  is  too  close  to  the  line  center.  With  the 
use  of  a  value  of  C  «  1 ,  it  could  become  necessary  to  renormalize  the  line 
strengths  in  order  to  handle  this  difficulty  consistently. 

A  quite  different  CO2  lineshape  model  has  actually  been  lmpl  -^ent^d  In 
FASC0D2.  This  was  done  partly  to  avoid  the  difficulty  of  fixing  the  "turn¬ 
on"  point  for  sub-Lorentzian  behavior,  and  partly  out  of  considerations 
related  to  the  computational  complexities  involved  in  modifying  all  of  the 
sub-functions,  The  sub-Lorentzian  form  of  the  lineshape  was  accounted  for 
by  multiplying  only  two  of  the  five  Lorentz  lineshape  subfunctions,  namely 
F4  and  the  REMAINDER  or  fifth  subfunction,  by  the  single  exponential-type 
form  factor  CHI  given  by 

B 

CHI (DELTA)  =  EXP(*  A  (DELTA)  ]  ,  [  4.5  ) 

and  A  and  B  are  adjustable  parameters.  Modification  of  only  the  last 
two  subfunctions  with  this  form  factor  produces  a  CO^  lineshape  given  by 

[  FI  +  F2  +  F3  ]  +  CHI  (DELTA)  [  F4  +  REMAINDER  ]  . 

From  Table  4.1,  it  is  clear  that  this  model  is  equivalent  to  a  lineshape 
specified  by 

L ( X )  -  [1  -  CHI (DELTA)]  Q64(X)  for  X  <  64 

and 

CHI (DELTA)  L(X)  for  X  >  64  . 
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This  lineshape  was  used  to  construct  a  number  of  trial  synthetic  CO2  ab¬ 
sorption  spectra  over  the  region  between  2360  and  2500  CM-1.  The  various 
continua  were  compared  to  laboratory  absorption  measurements  of  WSB  in 
order  to  adjust  A  and  B  for  the  best  fit.  The  values  A  =  0.623  and 
B  =  0.41  were  found  to  give  the  form  factor  which  is  closest  to  repro¬ 
ducing  the  laboratory  data  for  a  dilute  mixture  of  CO^  in  nitrogen  at 
T=296  K.  Table  4.2  lists  this  FASC0D2  form  factor  in  numerical  form  in 
addition  to  the  others  which  have  been  discussed  above. 

A  principal  advantage  of  the  FASC0D2  CO^  lineshape  is  that  the  problem 
of  specifying  a  turn-on  point  for  sub-Lorentzian  behavior  has  been  avoided. 

A  key  test  of  this  lineshape  is  whether  it  can  produce  a  continuum  absorp¬ 
tion  curve  which  compares  favorably  to  the  other  exponential -type  fits.  We 
have  found  the  FASC0D2  shape  to  fit  the  WSB  data  slightly  better  than  the 
best  simple  exponential  form  factor  with  C  *  0.5.  The  difference  is 
slight,  but  the  FASC0D2  lineshape  model  does  appear  to  be  more  than  ade¬ 
quate  to  fit  the  data.  A  first  comparison  was  made  with  the  simulated  C02 
continuum  transmittance  for  a  5.1  km  atmospheric  path  reported  in  [26]. 

The  FASC0D2  predictions  are  displayed  in  Fig.  4.1  as  points  on  the  trans- 
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mittance  curve  of  Roney  et  al,  The  comparison  ref  lects  a  difference  of 
only  1-2  percent  between  the  two  models. 

Further  comparisons  are  shown  in  Figs.  4.2  and  4.3.  In  Fig.  4.2  the 
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stratospheric  transmittance  data  of  Smith  et  al  is  reproduced  (a)  and 
displayed  above  a  FASC0D2  synthetic  spectrum  (b)  made  at  AFGL  using  the  new 
CO,,  continuum  model.  The  agreement  appears  to  be  quite  good,  with  FASC0D2 
most  nearly  matching  the  curve  for  the  lineshape  model  of  Susskind  and  Mo 
(curve  2)  at  the  bandhead.  Fig.  4.3  shows  the  6.4  km  (P  -  880  mbar,  T<*33°C) 
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path  transmittance  (a)  measured  by  Hanley  et  al  above  the  corresponding 
FASC0D2  calculation  (b).  Again  agreement  is  quite  good,  indicating  that 
the  FASC0D2  CK^  form  factor  does  give  an  improved  fit  to  the  far  wing 
absorption,  at  least  near  the  4.3  ym  bandhead. 


TABLE  4,2  COMPARISON  OF  VARIOUS  C02  LINESHAPE  FORM  FACTORS 


DELTA  I CM- 1 ] 

FASC0D2 

WSB 

BURCH 

SUSS. , MO 

CANN  ET  AL 

2.0 

0.895189 

1.000000 

0.414740 

0.488393 

0.811930 

4.0 

0.586513 

1.000000 

0.286610 

0.355176 

0.647770 

6.0 

0.283527 

0.631284 

0.208990 

0.283801 

0.507540 

8.0 

0.231923 

0.466504 

0.159790 

0.236701 

0.391230 

10.0 

0.201623 

0.381189 

0.127360 

0.202576 

0.298840 

12.0 

0.178058 

0.324358 

0.105270 

0.176458 

0.230360 

14.0 

0.159098 

0.282553 

0.089794 

0.155722 

0.185820 

16.0 

0.143463 

0.250046 

0.078712 

0.138817 

0.165190 

18.0 

0.130322 

0.223839 

0.070627 

0.124757 

0.156520 

20.0 

0.119111 

0.202165 

0.064643 

0.112875 

0.148310 

22.0 

0.109428 

0.183890 

0.060164 

0.102705 

0.140540 

24.0 

0.100979 

0.168248 

0.056780 

0.093905 

0.133190 

26.0 

0.093543 

0.154697 

0.054202 

0.086222 

0.126240 

28.0 

0.086949 

0.142838 

'.052221 

0.079462 

0.119660 

30.0 

0.081065 

0.132373 

0.050681 

0.073474 

0.113440 

32.0 

0.075783 

0.123070 

0.049462 

0.068138 

0.107570 

34.0 

0.071019 

0.114748 

0.048472 

0.063358 

0.102020 

36.0 

0.066702 

0.107264 

0.047638 

0.059056 

0.096774 

38.0 

0.062774 

0.100500 

0.046900 

0.055168 

0.091822 

40.0 

0.059188 

0.094360 

0.046213 

0.051640 

0.087148 

42.0 

0.055901 

0.088766 

0.045538 

0.048429 

0.082736 

44.0 

0  052881 

0.083650 

0.044848 

0.045495 

0.078572 

46.0 

0.050098 

0.078957 

0.044120 

0.042808 

0.074644 

48.0 

0.047526 

0.074639 

0.043340 

0.040340 

0.070939 

50.0 

0.045144 

0.070656 

0.042496 

0.038068 

0.067444 

52.0 

0.042933 

0.066972 

0.041585 

0.035971 

0.064148 

54.0 

0.040877 

0.063558 

0.040605 

0.034031 

0.061040 

56.0 

0.038960 

0.060387 

0.039558 

0.032233 

0.058110 

58.0 

0.037170 

0.057435 

0.038451 

0.030563 

0.055347 

60.0 

0.035496 

0.054683 

0.037290 

0.029011 

0.052743 

62.0 

0.033928 

0.052113 

0.036084 

0.027564 

0.050287 

64.0 

0.032457 

0.049709 

0.034844 

0.026213 

0.047972 

66.0 

0.031074 

0.047456 

0.033580 

0.024951 

0.045789 

68.0 

0.029774 

0.045342 

0.032302 

0.023770 

0.043730 

70.0 

0.028549 

0.043355 

0.031021 

0.022663 

0.041789 

72.0 

0.027394 

0.041487 

0.029746 

0.021625 

0.039958 

74.0 

0.026303 

0.039727 

0.028486 

0.020649 

0.038230 

76.0 

0.025271 

0.038067 

0.027248 

0.019731 

0.036599 

78.0 

0.024295 

0.036500 

0.026040 

0.018866 

0.035059 

80.0 

0.023371 

0.035019 

0.024866 

0.018052 

0.033605 

82.0 

0.022495 

0.033618 

0.023732 

0.017283 

0.032231 

84.0 

0.021663 

0.032292 

0.022640 

0.016557 

0.030933 

86.0 

0.020873 

0.031035 

0.021594 

0.015870 

0.029705 

88.0 

0.020122 

0.029842 

0.020595 

0.015221 

0.028543 

90.0 

0.019407 

0.028710 

0.019643 

0.014606 

0.027443 
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DELTA  [CM- 1 ] 

FASC0D2 

WSB 

92.0 

0.018727 

0.027634 

94.0 

0.018079 

0.026612 

96.0 

0.017461 

0.025638 

98.0 

0.016871 

0.024711 

100.0 

0.016308 

0.023828 

102.0 

0.015770 

0.022986 

104.0 

0.015255 

0.022182 

106.0 

0.014763 

0.021415 

108.0 

0.014292 

0.020682 

110.0 

0.013841 

0.019981 

112.0 

0.013400 

0.019311 

114.0 

0.012994 

0.018670 

116.0 

0.012596 

0.018056 

118.0 

0.012214 

0.017468 

120.0 

0.011848 

0.016905 

122.0 

0.011496 

0.016364 

124.0 

0.011158 

0.015846 

126.0 

0.010832 

0.015348 

128.0 

0.010519 

0.014871 

130.0 

0.010218 

0.014412 

132.0 

0.009928 

0.013971 

134.0 

0.009649 

0.013547 

136.0 

0.009380 

0.013139 

138.0 

0.009121 

0.012747 

140.0 

0.008871 

0.012370 

142.0 

0.008630 

0.012006 

144.0 

0.008397 

0.011657 

146.0  • 

0.008172 

0.011319 

148.0 

0.007956 

0.010994 

150.0 

0.007746 

0.010681 

152.0 

0.007544 

0.010379 

154.0 

0.007348 

0.010087 

156.0 

0.007159 

0.009806 

158.0 

0.006976 

0.009535 

160.0 

0.006800 

0.009273 

162.0 

0.006629 

0.009019 

164.0 

0.006463 

0.008775 

166.0 

0.006303 

0.008538 

168.0 

0.006147 

0.008310 

170.0 

0.005997 

0.008089 

172.0 

0.005851 

0.007875 

174.0 

0.005710 

0.007668 

176.0 

0,005573 

0.007468 

178.0 

0.005441 

0.007275 

180.0 

0.005312 

0.007087 

182.0 

0.005187 

0.006906 

184.0 

0.005066 

0.006730 

BURCH 

SUSS. ,M0 

CANN  ET  AL 

0.018739 

0.014023 

0.026401 

0.017883 

0.013470 

0.025413 

0.017074 

0.012945 

0.024477 

0.016310 

0.012446 

0.023587 

0.015591 

0.011972 

0.022742 

0.014913 

0.011521 

0.021939 

0.014276 

0.011092 

0.021175 

0.013677 

0.010683 

0.020447 

0.013114 

0.010294 

0.019754 

0.012584 

0.009922 

0.019092 

0.012087 

0.009568 

0.018460 

0.011619 

0.009229 

0.017857 

0.011179 

0.008906 

0.017280 

0.010765 

0.008597 

0.016728 

0.010374 

0.008302 

0.016199 

0.010005 

0.008019 

0.015692 

0.009656 

0.007749 

0.015206 

0.009326 

0.007490 

0.014739 

0.009012 

0.007241 

0.014291 

0.008714 

0.007003 

0.013861 

0.008431 

0.006775 

0.013447 

0.008160 

0.006557 

0.013049 

0.007901 

0.006347 

0.012665 

0.007653 

0.006145 

0.012297 

0.007414 

0.005951 

0.011941 

0.007185 

0.005765 

0.011599 

0.006964 

0.005587 

0.011270 

0.006750 

0.005415 

0.010953 

0.006543 

0.005249 

0.010647 

0.006343 

0.005090 

0.010353 

0.006149 

0.004937 

0.010070 

0.005960 

0.004790 

0.009797 

0.005776 

0 . 004648 

0.009535 

0.005598 

0.004511 

0.009283 

0.005424 

0.004379 

0.009040 

0.005255 

0.004252 

0.008807 

0.005090 

0.004129 

0.008584 

0.004929 

0.004011 

0.008370 

0.004773 

0.003897 

0.006164 

0.004620 

0.003787 

0.007967 

0.004472 

0.003681 

0.007779 

0.004329 

0.003578 

0.007599 

0.004189 

0.003479 

0.007428 

0.004053 

0.003383 

0.007265 

0.003922 

0.003291 

0.007109 

0.003795 

0.003201 

0.006961 

0,003672 

0.003115 

0.006821 
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DELTA  (CM- 1 ] 

FASC0D2 

186.0 

0.004948 

188.0 

0.004834 

190.0 

0.004723 

192.0 

0.004616 

194.0 

0.004511 

196.0 

0.004410 

198.0 

0.004311 

200.0 

0.004215 

202.0 

0.004122 

204.0 

0.004032 

206.0 

0.003943 

208.0 

0.003858 

210.0 

0.003774 

212.0 

0.003693 

214.0 

0.003614 

216.0 

0.003537 

218.0 

0.003463 

220.0 

0.003390 

222.0 

0.003319 

224.0 

0.003250 

226.0 

0.003183 

228.0 

0.003117 

230.0 

0.003053 

232.0 

0.002991 

234.0 

0.002930 

236.0 

0.002871 

238.0 

0.002813 

240.0 

0.002757 

242.0 

0.002702 

244.0 

0.002649 

246.0 

0.002597 

248.0 

0.002546 

250.0 

0.002496 

WSB  BURCH 

0.006560  0.003553 
0.006395  0.003439 
0.006235  0.003329 
0.006080  0.003223 
0.005930  0.003121 
0.005784  0.003024 
0.005642  0.002930 
0.005505  0.002841 
0.005372  0.002755 
0.005243  0.002674 
0.005118  0.002596 
0.004996  0.002522 
0.004878  0.002452 
0.004763  0.002385 
0.004651  0.002322 
0.004543  0.002262 
0.004437  0.002205 
0.004335  0.002151 
0.004236  0.002099 
0.004139  0.002051 
0.004045  0.002005 
0.003953  0.001961 
0.003864  0.001919 
0.003778  0.001880 
0.003693  0.001842 
0.003611  0.001806 
0.003532  0.001771 
0.003454  0.001738 
0.003378  0.001706 
0,003305  0.001675 
0.003233  0.001645 
0.003163  0.001615 
0.003095  0.001586 
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0.003031  0.006689 
0.002950  0.006563 
0.002872  0.006445 
0.002796  0.006334 
0.002723  0.006229 
0.002652  0.006131 
0.002583  0.006040 
0.002517  0.005954 
0.002452  0.005875 
0.002390  0.005801 
0.002329  0.005733 
0.002270  0.005670 
0.002213  0.005613 
0.002158  0.005560 
0.002105  0.005513 
0.002053  0.005469 
0.002002  0.005430 
0.001953  0.005395 
0.001906  0.005364 
0.001860  0.005336 
0.001815  0.005311 
0.001772  0.005290 
0.001729  0.005271 
0.001688  0.005255 
0.001649  0.005241 
0.001610  0.005229 
0.001572  0.005219 
0,001536  0.005210 
0.001500  0.005203 
0.001466  0.005197 
0.001432  0.005192 
0.001399  0.005188 
0.001367  0.005184 


Figure  4.3  Measured  atmospheric  transmittance  over  a  6.4  Km  path  for  8.8  torr  H^O  and 
33°C  air  temperature  as  reported  by  Hanley  et  al^7  (a)  as  compared  to  (b) 
a  corresponding  FASC0D2  spectrum. 


5.  LAYERING  IN  FASCOD 


Radiative  transfer  in  FASCOD  is  calculated  in  a  special  medium  com¬ 
posed  of  a  number  of  homogeneous  layers  that  approximate  the  inhomogeneous 
atmosphere.  The  manner  in  which  these  layers  are  choosen  and  how  their 
radiative  properties  are  obtained  effect  the  accuracy  of  the  radiative 
transfer.  This  section  discusses  this  problem  by  first  developing  the 
Curtis-Godson  approximation  which  is  a  two-parameter,  equivalent- line 
argument  for  inhomogeneous  media.  However,  such  an  approximation  is 
cumberson  to  execute  in  the  context  of  FASCOD.  An  approximation  to 
Curtis-Godson  Is  therefore  developed  that  leada  to  what  is  presently 
done  in  the  various  versions  of  FASCOD.  With  this  information  available, 
the  last  section  makes  recommendations  for  further  work. 

5. 1.  Curtis-Godson  Approximation:  Radiative  transfer  through  in¬ 
homogeneous  media  is  always  difficult  to  model  in  that  all  the  associated 
spatial  integrals  must  be  evaluated  numerically  on  a  fine  grid.  This  is 
operationally  cumbersome  and  computationally  expensive.  Consequently 
several  approaches  have  been  developed  for  replacing  inhomogeneous  media  by 
homogeneous  media  that  are  in  some  way  radlatlvely  equivalent.  The  basis 
for  most  of  these  methods  is  the  Curtis-Godson  approximation  which  is  a 
two-parameter  equivalent  line  argument.  This  approximation  may  be  developed 
from  several  different  viewpoints,  but  the  common  principle  is  that  certain 
radiative  transfer  properties  are  matched  to  find  the  required  parameters. 

The  exact,  single-line,  spectral  optical  depth  along  some  path 
segment  in  a  given  layer  for  a  specific  absorber  n  can  be  expressed  as 

/Z2  du 

Sn(T)fvn(T,P)  (~)  dz  .  (5.1) 

Zi 
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Each  layer  contains  N  different  molecular  absorbers  designated  by  a 
counting  index  n  with  total  column  densities  or  amounts  un  per  unit 
cross-sectional  area  given  by 


(5.2) 


where  Z  is  the  vertical  dimension.  The  derivative  includes  the  density 
times  the  secant  of  the  local  zenith  angle  (that  may  vary  with  Z)  so  the 
path  length  is  accounted  for.  The  line  shape  function  f  is  normalized 
such  that 


f(v)dv 


1. 


(5.3) 


and  S  is  the  line  strength.  An  equivalent  homogeneous  medium  would  have 
an  optical  depth  given  by 


where  T^n  is  again  normalized  as  in  Eq.  (5.3).  The  spectrally  averaged 
optical  depth  of  expressions  (5.1)  and  (5.4)  are  equated  (which  is  the 
same  as  saying  the  two  approaches  are  equal  in  the  weak-line  limit)  which 
gives  the  relation 


du 

Sn(T)  <—)  dZ  . 


(5.5) 


pi 


47 


This  says  that  the  line  strength  to  use  in  the  homogeneous  layer  is  S 
If  the  radiative  transfer  is  also  forced  to  match  in  the  strong-line 
limit  where  the  line  wings  dominate,  then  fv  takes  the  form 


a 


v  <v-V‘ 


(5.6) 


where  a  is  the  line  halfwidth.  One  therefore  obtains  the  equality  that 


*  n 
a  u 

c  n  _ n 

&  2 

(v-v0)2 


(v-ve)- 


I  s"tt>  “  “jr 


dZ 


(5.7) 


or  the  equivalent  halfwidth  is  of  the  form 

*1, 


5° 


f  2  /du  \ 

l  Sn(T)  a(T,  P)^j 


dZ 


j2  s"«)(h^)« 


(5.8) 


-1/2 

Since  a  is  proportional  to  P(T)  ,  one  obtains 

•Z„  y m  v 1/2, 


5n 


a 


l  •  m 


sn 


(5.9) 


where  a  is  a  reference  halfwidth  and  Pn  and  T  are  the  reference  pressure 
o  o  o 

and  temperature,  respectively. 
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The  Curtis-'Godson  approximation  replaces  any  given  molecular  line  in 
an  inhomogeneous  medium  with  a  "radiatively  equivalent"  line  in  an  homogeneous 
medium  that  has  a  single  linestrength  ?  n  and  halfwidth  a  n.  The  accuracy 
of  this  approximation  is  not  easily  qualified  since  it  does  depend  on  the 
specified  inhomogeneous  structure.  To  check  the  accuracy,  numerical  experi¬ 
ments  are  usually  performed. 

The  Curtis-Godson  approximation  says  that  strength-weighted  averages 
should  be  used.  This  would  be  a  complicated  and  time  consuming  approach 
since  such  averages  would  have  to  be  performed  for  each  line.  The  next 
section  develops  a  further  approximation  of  the  Curtis-Godson  approach 
that  is  much  more  efficient  to  execute  in  FASCOD. 

5.2  Present  Layering  Procedure:  If  the  line  strength  Sn(T)  in  any  layer 
between  and  in  the  atmosphere  is  expanded  in  a  Taylor  series  about 
the  amount  weighted  temperature  T  given  by 

dZ  ,  (5.10) 

and  Eq.  (5.5)  gives 

;  [2(¥)2(~si)dz  *  (5-u> 

'  "h 

FASCOD  uses  Sn(T)  as  the  layer's  line  strength  and  this  agrees  with  the 
Curtis-Godson  approximation  to  the  second-order  term  shown.  For  uniformly 
mixed  gases,  one  has  that 


t  /'"“’(S') 


then  the  linear  term  is  identically  zero 

l 


n.2_n 

s  »  .  S"(T)  + 

8i 


_1 
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(5.12) 


where 


and 


-  C 


-£ 


u 

n 


1 


(5.13) 


(5.14) 


Therefore  the  temperature  T  can  also  be  expressed  as 


T  - 


> 


(5.15) 


which  is  the  mean  temperature  actually  used  In  FASCOD.  The  reader  shc-’d  be 
aw Are,  however,  that  FASCOD  includes  nonuniform  gases  in  the  sum  of  Eq. 
(5.13)  and  this  is  not  consistent  with  the  present  development.  To 
complete  the  argument,  Sn(T)  and  a(T,P)  are  expanded  in  Taylor  series 
in  Eq.  (5.8)  about  T  and  P,  where  P  is  given  by 

Z2 

%  /  P(z) 

Z1 

Again  the  linear  tenus  will  be  identically  zero  and  one  obtaines 


dZ 


(5.16) 


—  n 
a 


an(T,  P)  +  second-order  terms  , 


(5.17) 
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where  again  the  mean  pressure  P  used  in  FASCOD  is 


(5.18) 


Using  the  imperical  relationship  stated  above  Eq.  (5.9)  for  a  leads  to  the 
expression 


To\  54 

—  )  +  second-order  terms  (5.19) 

T  / 

that  is  used  in  FASCOD.  The  conclusion  is  that  FASCOD' s  layer  properties 
are  obtained  by  a  Curtis-Godson  approximation  valid  to  second  order  in 
terms  like  (T-T)2  or  (P-P)2  [see  Eq.  (5.11)].  The  use  of  S(T)  and  a(T,  P) 
for  each  layer  in  the  atmosphere  is  computationally  very  efficient  since 
line-by-line  averages  need  not  be  performed. 

Once  the  radiative  properties  for  each  layer  are  properly  obtained, 
the  problem  of  choosing  the  layer  boundaries  must  be  faced.  All  the 
versions  of  FASCOD  used  in  the  work  discussed  in  this  report  have  the 
layer  boundaries  specified  by  the  user.  The  only  requirement  is  that  the 
average  line  halfwidths  cannot  vary  by  more  than  the  ratio  of  2  to  1  in 
adjacent  layers.  A  check  on  this  condition  is  performed  and  execution 
terminated  when  the  requirement  is  not  met.  This  is  inconvenient,  and 
later  versions  of  FASCOD  now  active  at  AFGL  create  more  layers  and  adjust 
boundaries  automatically  until  the  above  criterion  (or  any  other  value  of 
this  ratio  set  by  the  user)  is  satisfied.  These  later  versions  of  FASCOD 
also  contain  a  restriction  on  the  temperature  difference  allowed  for  adjacent 


51 


boundaries.  Both  of  these  constraints  keep  the  second-order  error  terms 
small  in  the  Curtis-Godson  arguments  given  above.  They  also  play  a  role  in 
maintaining  spectral  accuracy.  FASCOD  has  been  run  for  both  32  and  16 
layers  for  a  few  specific  problems  and  the  transmittance  and  radiance 
have  shown  very  little  change.  This  says  the  layering  procedure  is  very 
conservative  and  accurate. 

5.3.  Recommendations  For  Further  Studies:  As  brought  out  in  the  above 
discussion,  the  layering  in  FASCOD  is  conservative,  leads  to  accurate  results, 
and  the  radiative  properties  for  each  layer  are  being  calculated  properly. 

The  one  obvious  and  minor  change  to  be  made  is  that  the  nommlformly  mixed 
gases  should  be  handled  separately  in  defining  their  mee’*  temperature  and 
pressure.  Some  test  runs  could  be  performed  to  see  whether  this  is  an 
important  point  in  practice. 

Since  the  execution  time  of  FASCOD  is  dependent  on  the  number  of  layers 
used,  one  would  like  to  know  the  relationship  between  the  number  of  layers 
used  and  the  accuracy  of  the  results.  This  is  not  a  simple  question  due  to 
the  fact  that  the  code  operates  with  a  large  number  of  parameters,  all  of 
which  can  effect  the  answer.  As  mentioned  above, some  experiments  have  been 
done  on  changing  the  number  of  layers  with  a  fixed  set  of  parameters. 

Further  numerical  experiments  should  be  carried  out  so  the  layering  procedure, 
accuracy,  and  parametric  space  can  be  related  to  develop  practical  guidelines 
for  the  users  of  FASCOD. 
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6.  CONCLUDING  REMARKS 


As  observed  in  reading  this  report,  the  generic  code  called  FASCOD 
has  been  under  continuous  development  by  AFGL  and  Sonicraft  during  the 
period  covered  by  this  contract.  The  FASCOD  versions  used  varied 
throughout  the  work  and  this  was  noted  in  the  report.  Only  those  parti¬ 
cular  sections  involved  in  the  present  modifications  are  documented  here. 

The  developments  presented  here  plus  those  accomplished  at  AFGL 
will  be  combined  in  the  near  future  into  one  integrated  code  with  Increased 
capabilities.  This  integrated  code  will  be  released  to  the  public  as 
FASC0D2. 
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APPENDIX  A:  FASCOD2-80  USER'S  GUIDE 

This  appendix  documents  only  those  parts  of  the  code  that  have  been 
changed  to  include  the  new  geometry  routine,  the  gaseous  and  aerosol 
atmospheric  models,  and  the  integration  and  automation  of  several  programs. 
An  example  run  is  also  presented  and  the  output  format  explained. 

A.l  Input  Card  Sequence  and  Format;  In  the  course  of  adding  the  geometry 
and  aerosol  routines  to  FASCODE,  several  new  control  cards  were  Introduced. 
Some  of  the  original  FASC0D1A  control  cards  (same  as  FASC0D1B)  have  been 
consolidated  to  simplify  the  input  deck  structure.  The  complete  input 
card  sequence  is  listed  below.  Only  the  new  input  parameters  are  defined 
here.  The  user  is  referred  to  Ref.  [6]  for  a  description  of  the  original 
FASCODl  input  variables. 

Card  1  from  FASCODl  main 

XID(I),  I  -  1,  7 
(10A10) 

Card  2  from  FASCODl  main 

IHIRAC,  ILBLF4 ,  ICNTNM,  IAERSL,  I EMIT,  I SCAN,  IFILTR ,  IPLOT, 
IATM,  MPTS,  NPTS  (9(4X,  II),  25X,  215) 


where 


IAERSL  -  0 

-  1 

IATM  -  0 

-  1 


aerosol  attenuation  not  included 

aerosol  attenuation  included 

ATMPTH  geometry  routines  are  used 

ATMPTH  is  not  used;  layer  data  is  read  from  TAPE7 
as  before 


Card  3  from  FASCODl  main 


VI,  V2,  TBOUND ,  EMISIV  (10E10.3) 


Card  4  from  AERSOL,  included  only  if  aerosol  attenuation  is  to  be  included, 
i.e.,  IAERSL  -  1 

I HAZE,  ISEASN,  1VULCN,  JP,  VIS  (4  (4X,  II),  5X,  E10.3) 

where, 

IHAZE  -  0  no  aerosol  attenuation  included  in  the  calculation. 

-  1  RURAL  extinction,  23  km  VIS. 

■  2  RURAL  extinction,  5  km  VIS. 

-  3  MARITIME  extinction,  23  km  VIS. 

-  4  MARITIME  extinction,  5  km  VIS. 

-  5  URBAN  extinction,  5  km  VIS. 

-  6  TROPOSPHERIC  extinction,  50  km  VIS. 

“  7  USER  DEFINED  extinction,  23  km  VIS.  (Read  into  the  program 

immediately  after  CARD1.  Refer  to  the  main  program  LOWEM 
in  Appendix  A5  for  the  input  format  of  coefficients) . 

-  8  F0G1  (advection  fog)  extinction,  0.2  km  VIS. 

-  9  FOG 2  (radiation  fog)  extinction,  0.5  km  VIS. 

ISEASN  -  0  season  determined  by  the  value  of  MODEL; 

SPRING-SUMMER  for  MODEL  -0,  1,  2,  4,  6,  7 
FALL-WINTER  for  MODEL  -  3,  5 

-  1  SPRING-SUMMER 

-  2  FALL-WINTER 

IVULCN  -  0,  1  BACKGROUND  TROPOSPHERIC  profile  and  extinction 

-  2  MODERATE  VOLCANIC  profile  and  AGED  VOLCANIC  extinction 

■  3  HIGH  VOLCANIC  profile  and  FRESH  VOLCANIC  extinction 

-  4  HIGH  VOLCANIC  profile  and  AGED  VOLCANIC  extinction 

-  5  MODERATE  VOLCANIC  profile  and  FRESH  VOLCANIC  extinction 

JP  «  print  option  parameter,  inactive  in  this  version  of  the  code 


VIS  -  meteorological  range  (km) 

(when  specified,  supercedes  default  value  set  by  IHAZE) . 
Cards  5,  6,  7,  and  the  layer  boundary  cards  pertain  to  ATMPTH, 
and  IATM  •  0 

Card  5  MODEL,  ITYPE,  IIN,  IMOD,  KMAX,  RE  (515,  5X,  F10.4) 


where , 


MODEL  -  0: 

1: 

2: 

3: 

4. 

5: 

6: 

7s 


USER  SUPPLIED  horizontal  path  parameters 

TROPICAL  model  atmosphere 

MIDLATITUDE  SUMMER 

MIDLATITUDE  WINTER 

SUBARCTIC  SUMMER 

SUBARCTIC  WINTER 

U.S.  STANDARD,  1962 

USER  SUPPLIED  atmospheric  profile 


ITYPE  -  1:  horizontal  path  (constant  pressure) 
2:  slant  path  from  HI  to  H2 
3:  slant  path  from  HI  to  space  (100  km) 


IIN  ■  number  of  boundary  altitudes  for  the  FASC0D1  layers  ( required 
for  ITYPE  -  2  OR  3) . 

IMOD  ■  number  of  boundary  altitudes  for  USER  SUPPLIED  atmospheric 
profile  (MODEL  7),  default  ■  34  . 

KMAX  -  number  of  molecular  species  for  which  the  amounts  are  to 
be  calculated.,  default  ■  8  . 

RE  -  radius  of  the  earth,  defaults:  MODEL  «  1,  RE  -  6378.39  km; 
MODELS  -  2,  3,  6,  7,  RE  -  6371 . 23km;  MODELS  "  4,  5,  RE  - 
6356.91  km. 


58 


The  formats  of  control  cards  6  and  7  depend  on  whether  the  path  is 
horizontal  ( ITYPE  -  1)  or  slant  (ITYPE  -  2  or  3)  . 

For  a  slant  path: 

CARD  6:  Hi,  H2,  ANGLE,  RANGE,  BETA,  LEN 
(5F10.4,  15) 

HI  altitude  of  the  observer  or  receiver  (km) 

H2  altitude  of  the  other  endpoint  of  the  path  (km) 

ANGLE  zenith  angle  at  HI  to  H2  (km) 

RANGE  length  of  the  path  from  HI  to  H2 

BETA  earth  centered  angle  for  the  path  Hi  toH2  (degrees) 

LEN  -  0,  short  path;  1,  long  path  through  a  tangent  height. 

LEN  is  used  only  when  ANGLE  is  GT  90.0  and  HI  >  H2. 
default  -  0 

Only  three  of  the  first  five  parameters  need  be  specified;  for  example, 

HI.  H2,  ANGLE,  or  HI,  H2,  BETA,  or  HI,  ANGLE,  RANGE.  Refer  to  the  comments 
in  subroutine  GMTRY  or  see  Reference  [5]  for  a  discussion  of  allowable  com¬ 
binations  of  these  parameters.  Next  the  boundary  altitudes  for  the  FASCODE 
layers  are  read  in  with  format  (8F10.3). 

Card  7:  VI,  V2  (2F10.3) 

VI,  V2  initial  and  final  wavenumbers  for  use  in  calculating  the 
index  of  refraction  (cm-1) 

IF  MODEL  «  7,  the  input  atmospheric  profile  is  read  in  after  control  card  1 
in  the  .following  format 

(HEADER  (I),  I  -  1,  2)  (2A10) 

A  20  character  header  describing  the  profile  Z,  P,  T  ((3F10.3) 
(DENSITY  (K),  K  -  1,  KMAX)  (8E10.3): 


two  cards  for  each  of  the  IMOD  levels  giving  the  altitude 
(km),  pressure  (mb),  temperature  (K) ,  and  densities  of  the 
molecular  species  (molecular  cm-3)  at  each  level) 


For  a  horizontal  path; 

CARD  2:  Z,  P,  T,  RANGE  (DEN(K),  K  -  1,  KMAX)  (4F10.3,/,  (8E10.3) 

Z  altitude  (km) 

P  pressure  (mb) 

T  temperature  (K) 

RANGE  path  length  (km) 

DEN(K)  density  of  the  K'th  molecular  species  (molecules  cm-3) 

For  MODEL  ■  1  to  7,  only  Z  and  RANGE  are  used  and  P,  T,  and  DEN  are  interpolated. 
CARD  3:  NOT  USED 


For  model  7,  the  input  model  atmosphere  is  read  in  after  control  card  5  as  for 
a  slant  path. 

For  model  7  horizontal  or  slant  paths  the  aerosol  amounts  must  be  supplied 
if  aerosol  attenuation  is  to  be  included,  i.e.  IAERSL  ■  1.  The  following 
card  is  repeated  once  for  each  atmospheric  layer 

AWKAER(J)  FORMAT (E10. 3) 

A. 2  Sample  Run  of  FASCOD2-8Q 

The  sample  run  of  FASCOD2-80  using  the  geometry  and  aerosol  options 
presented  here  will  illustrate  the  complete  input  deck  structure  and  provide 
"standard"  output  for  comparison  with  the  user's  version  of  the  code.  The 
complete  input  file  follows. 
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SCRAF,  CM70000,  T60 


2710  NAME 


ATTACH,  TAPE3,  LINES2050T02150,  ID-NAME,  MR-1. 
ATTACH,  LGO,  FASCOD1NAME,  ID-NAME,  MR-1. 
REQUEST,  TAPE6,  *PF. 

MAF,  PART. 

LDSET,  PRESET-INDEF. 

LOAD,  LGO. 

SEGLOAD. 

EXECUTE,  FASCOD1. 

EXIT(U) 

CATALOG  (TAPE6,  FASOUT,  ID-NAME) 


EOR 


FASCODE  TREE  FASCOD1  -  (PATH,  HIRAC1,  LBLF4,  CONTNM,  ABSMRG,  EMINIT,  RADMRG, 
,  SCANFN ,  PLOTT,  AERSOL,  ATMPTH) 

FASCOD1  GLOBAL  MAIN,  CONSTS,  L1NHDR ,  FILHDR,  ABSORB,  SCATTR ,  LBLF4,  PROF, 

,  MISC,  OUTPUT 

LBLF4  EQUAL  RFRPTH 

PATH  INCLUDE  PATH 

PATH  GLOBAL  PATH-SAVE 

HIRAC1  INCLUDE  HIRAC1,  SHAPEL,  SHAPEG,  VO ICON,  RDL1N,  CNVFNV,  PANEL,  MOLEC, 

,  QV,  CONT,  ABSOUT 

LBLF4  INCLUDE  LBLF4 ,  RDLIN4,  CONVF4 ,  MOLEC,  QV,  SHRINK 


CONTNM 

INCLUDE 

CONTNM, 

XINT 

ABSMRG 

INCLUDE 

ABSMRG, 

ABSOUT 

EMINIT 

INCLUDE 

EMINIT, 

EMIN,  EMOUT,  BBFN 

RADMRG 

INCLUDE 

RADMRG, 

EMIN,  EMCUT ,  BBFN 

SCANFN 

INCLUDE 

SCANFN, 

SHAPEG 

AERSOL 

INCLUDE 

AERSOL 

AERSOL 

GLOBAL 

ICALL , 

COEFF,  AER-SAVE 

PLOTT  TREE  PLTFAS  -  (HEADER,  AXEST,  FPLINE) 

AXEST  TREE  AXES  -  (AXISL,  AXLOG,  AX2) 

END  FASCOD1 

DATA  FOR  FASCOD2  REPORT,  GROUND  TO  20  KM,  VERTICAL  AEROSOLS  ON 


HI  -  1 

F4 

-  1  CN  -  1 

AE  -  1 

EM  -  1  SC  -  1 

FI 

-  0 

PL  -  0 

2.095E  +  03 

2.105E  +  03 

0.000E  +  00  0.000E  +  00 

0.000E  +  00 

1H  -  1 

IS 

-  1  IV  -  1 

JP  -  1 

2.30E  +  01 

6 

2 

21 

00.0 

20.0 

00.0 

0. 

1. 

2. 

3.  4. 

5. 

6. 

7. 

8. 

9. 

10. 

11.  12. 

13. 

14. 

15. 

16. 

17. 

18. 

19.  20. 

21. 

22. 

23. 

2090. 

2110. 

1. 

2095.1 

2104 

.95  0 

3 

-1. 

EOR 

EOF 

This  file  consists  of  three  logical  records.  The  first  record  contains 
the  job  control  language  necassary  to  accesB,  load  and  execute  the  program. 
The  second  record  contains  the  directives  used  to  segment  the  program 
during  loading  to  minimize  the  central  memory  requirements.  The  third 
record  contains  the  input  data  described  earlier.  Note,  the  precise  spacing 
of  data  Is  not  preserved  in  the  sample  input  file  shown  above.  The  user 
is  referred  back  to  Section  A.l  for  the  exact  format  of  the  data.  This 
particular  run  la  in  the  emission  mode  using  the  new  geometry  and  aerosol 
routines  (IATM  ■  0  and  IAERSL  "1).  A  vertical  path  from  0  to  20  km  was 
chosen,  with  layer  boundaries  every  kilometer.  The  scanning  routine  is 
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utilized  (ISCAN  ■  1)  while  the  plotting  routine  which  is  machine  de¬ 
pendent  is  not  (IPLOT  ■  0) . 

The  following  output  file  has  been  truncated  to  include  only  output 
up  to  the  second  layer,  plus  the  final  layer  and  the  results  of  the 
scanning  routine.  This  should  be  sufficient  for  a  user  to  verify  that 
his  copy  of  the  code  is  functioning  properly.  The  basic  format  of  the 
output  file  remains  unchanged.  Additional  output  Includes  the  following: 

1.  A  summary  of  the  AERSOL  control  parameters  as  the  fourth  line 
of  output. 

2.  A  list  of  the  altitude  dependent  aerosol  densities,  labelled 
as  "AERSOL  PROFILE"  on  the  outfile. 

3.  A  table  of  aerosol  extinction  and  absorption  coefficients  In 
each  of  the  four  altitude  regions  as  a  function  of  wavelength. 
This  table  is  labelled  "EXTINCTION  AND  ABSORPTION  COEFFICIENTS", 
and  appears  about  midway  through  the  output  for  layer  number  one. 

4.  The  Integrated  aerosol  amount  for  each  layer,  labelled 
"WKAER-". 


DATA  r on-  PASC0P2  REPORT.  GROUND  TO  20  KMr  VERTICAL r AEROSOLS  ON  “A5C0D1  80/03/01.  14.43.04 
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APPENDIX  B:  FASC0D2-82  PROGRAM  REVISIONS  AND  USERS  GUIDE 


The  code  designated  FASCOD2-82  is  based  on  the  AFGL  code  FASCOD1B  of 
May,  1981.  Modifications  have  been  made  in  order  to  (1)  adjust  the  far 
wing  C0„  lineshape  using  a  new  fourth  subfunction  form  factor  CHI  and  a  new 
continuum  and  (2)  allow  non-LTE  line  strengths  to  be  used  at  high  altitude 
in  both  transmittance  and  radiance  calculations. 

In  the  case  of  CO2,  the  fourth  function  lineshape  is  now  computed 
separately  from  the  other  gases.  The  same  program  lineshape  control  para¬ 
meters  apply.  These  govern  which  of  the  Lorentz  subfunctions  are  used,  and 
the  various  combinations  are  enumerated  below: 

IHIRAC=1  allows  the  line-by-line  summation  algorithm  to  be  used 
with  the  three  subfunctions  FI,  F2,  and  F3  which  cover  the  region 
within  64  halfwidths  of  line  center  (for  all  molecules). 

ILBLF4»1  adds  the  fourth  subfunction  F4  (if  IHIRAC*®1)  to  the 
first  three  subfunctions.  F4  is  modified  by  the  exponential 
form  factor  CHI  in  the  case  of  CO2  only. 

ICNTNM=1  adds  the  scaled  "fifth"  or  remainder  continuum  sub¬ 
function  for  CO2,  N9 ,  and  water  vapor. 

The  following  summary  shows  the  structure  of  program  routines  related  to 
the  fourth  and  "fifth"  subfunctions.  The  asterisk  indicates  where 
FASC0D2-82  routines  contain  significant  changes.  These  program  changes  for 
the  sub-Lorentzian  C07  lineshape  have  been  largely  confined  to  the  routines 
SHRINK,  C0NVF4,  and  the  CO^  continuum  data  subroutine  FRNC02. 
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LBLF4  ROUTINES 


LBLF4  Fourth  function  driver.  Scales  line  strengths  by 

column  density  amounts  and  temperature  correction 
factors.  Calculates  P  and  T  adjusted  halfwidths 
ALFAD  and  ALFAV.  Modifies  R4  by  radiation  field 
term  and  returns  R4  to  LAYER  for  later  merger  with 
ABSRB,  R3 ,  R2 ,  and  R1 . 

VOICON  Computes  the  Voigt  profile  factor  AVRAT(IZETA)  . 

MOLEC  Computes  line  strength  (SCOR)  and  Lorentz  halfwidth 

(ALFCOR)  correction  factors  from  P  and  T.  Also  cal¬ 
culates  the  Doppler  halfwidth  (ALFD1)  at  GNU=1. 

RDLIN4  Up  to  1250  lines  are  loaded  into  the  arrays  GNU(I), 

S(I),  ALFAL(I),  ALFAD(I)  and  MOL(I). 

*  The  Non-LTE  vibrational  transition  index  is  stripped 
from  the  MOL  array. 

BUFIN  is  called  to  read  blocked  line  data  from  TAPE3. 

SHRINK  Merges  nearby  lines  together  to  form  fewer  effective 

lines . 

*  CO2  and  noti-C02  lines  are  'lumped'  separately. 

Returns  the  new  line  data: 

GNU(J)  strength-weighted  average  frequency 

S(J)  summed  line  strength  for  merged  lines 

ALFAL(J)  strength-weighted  Lorentz  haltwidth 

ALFAD(J)  strength-weighted  Doppler  halfwidth 

*  MOL(J)  mo. ocular  ID  “0  for  non-C02  lines 

=2  for  C02  lines 

C0NVF4  Performs  line-by-line  fourth  function  sum  to  rroduce 

the  R4  array. 

Uses  a  Lorentz  subfunction  F4  which  is  modified  by 
^he  exponential  form  factor  CHI  for  CO^  lines  only. 

RACFN  Returns  the  value  of  the  radiation  field  term. 


CONTINUUM  ROUTINES 


CONTNM 

Scales/ interpolates  stored  continuum  data  and  loads 
the  array  ABSRB. 

SLF296 

Stored  self-broadened  H  0  continuum  for  T=296  K. 

2 

SLF260 

Stored  self-broadened  H  0  continuum  for  T=260  K. 

2 

FRN296 

Stored  foreign-bioadened  H^O  continuum  for  T=296  K. 

FRNC02 

>v  Stored  foreign-broadened  CO  2  continuum  for  T=296  K. 
Includes  the  sub-Lorentz  'REMAINDER'  subfunction. 

CONT 

Called  by  HIRAC  to  merge  ABSRB  and  R4  into  R3. 

The  non-LTE  calculation  uses  a  straightforward  extension  of  the  HIRAC 
line-by-line  algorithm,  which  is  made  to  perform  two  line-by-line  convolu¬ 
tions  simultaneously.  An  entirely  new  program  branch  (bearing  the  interim 
name  HIRACQ)  is  activated  at  high  altitudes  to  perform  the  non-LTE  calcula¬ 
tions  if  the  user  specifies  INLTE*!.  The  new  branch  routines  HIRACQ, 

CNVFNQ ,  and  PANELQ  are  the  non-LTE  analogs  of  HIRAC1,  CNVFNV,  and  PANEL. 

The  new  subroutine  STRTHS  is  used  by  HIRACQ  to  compute  the  two  "effective 
strengths"  for  each  line  which  is  enhanced  or  diminished  by  the  non-thermal 
vibrational  populations.  The  vibrational  populations  of  H2O,  CO2,  O2,  and 
NO  are  functions  of  altitude  and  are  read  directly  from  the  appropriate 
output  file  of  the  Degges  and  Smith  high  altitude  computer  model.  The 
populations  read  from  this  file  are  used  to  calculate  the  vibrational  state 
population  enhancement  ratios  by  the  subroutine  RDPOPS.  The  ratios  are 
stored  for  use  by  STRTHS.  The  actual  single-layer  radiance  calculation  is 
performed  by  EMIN,  which  has  been  modified  to  calculate  radiance  from  the 
two  HIRACQ  line-by-line  absorption  and  "effective  absorption"  spectral 
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functions.  A  summary  of  these  routines  with  brief  descriptions  are  given 
below: 


NON-LTE  ROUTINES 


BCDNLTE 

RDPOPS 

HIRACQ 

STRTHS 

CONVFNQ 
PANE  I, Q 


EMIN 


Prepares  a  blocked  line  file  for  use  in  the  non-LTE 
calculation.  Vibrational  state  labels  are  coded  and 
added  to  the  molecular  ID  for  each  non 'equilibrium  line. 

Reads  the  file  containing  non-LTE  population  data. 
Calculates  population  enhancement  ratios. 

Performs  two  separate  line-by-line  calculations. 

Computes  the  effective  line  strength  for  absorption 
and  for  the  second  function  C.  Returns  these  two 
strengths  to  HIRACQ  for  each  line. 

Constructs  the  arrays  (R1,R2,R3;  RR1,RR2,RR3)  which 
are  needed  for  the  two  line-by-line  sums. 

Merges  the  two  separate  sets  of  subfunctions  and 
writes  the  spectral  functions  K  and  C  to  a  temporary 
file. 

Calculates  the  single  layer  radiance  and  transmittance 
from  the  spectral  functions  K  and  C. 


The  off-line  program  BCDNLTE  prepares  the  blocked  line  data  file  for 
use  by  FASC0D2  in  doing  non-LTE  calculations.  In  addition  to  formatting 
the  line  data  (frequency,  strength,  halfwidth,  lower  state  energy,  and 
molecular  ID),  this  program  sorts  the  linos  into  equilibirum  and  non-equil¬ 
ibrium  transitions,  and  then  adds  two  vibrational  state  code  numbers  to  the 
molecular  ID.  The  code  numbers  are  presently  based  on  the  order  of  states 
used  in  the  model  of  Degges  and  Smith,  but  can  be  altered  to  allow  use  with 
other  non-LTE  gas  models  in  the  future.  More  details  about  program  BCDNLTE 
can  he  found  in  Appendix  C. 
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On  the  next  two  pages  a  sample  job  submission  file  is  presented.  This 
file  can  be  used  to  compile  and  execute  FASCOD2-82  if  the  molecular  data 
has  been  prepared  previously  by  BCDNLTE.  Output  from  BCDNLTE  has  been 
stored  in  file  "NLTELINEFILE"  in  this  example,  File  "BASEVIRPOP" ,  the 
Degges  and  Smith  output  file,  contains  the  needed  vibrational  population 
data. 


The  listing  below  illustrates  two  job  features.  The  SEGLOAD  tree  can 
be  used  to  segment  the  program  for  execution  with  about  70K  of  space  on  the 
AFGL  CDC  6600  computer.  The  tree  shows  the  heirarchy  of  program  subrou¬ 
tines  and  the  disposition  of  various  COMMON  blocks.  The  last  eight  lines 
also  present  a  sample  job  INPUT  card  sequence.  The  program  control  param¬ 
eters  and  data  which  appear  will  not  be  described  in  detail  here  since  only 
one  change  has  been  made  to  the  corresponding  card  sequence  for  version 
FASC0D1B.  The  single  additional  parameter  (INLTE)  appears  at  the  end  of 
card  2  (EQ=0).  In  this  case  the  non-equilibrium  option  has  been  turned 
off,  so  that  HIRACQ  will  not  be  executed  since  the  run  is  to  be  made  for  a 
(equilibrium)  sea  level  horizontal  path. 


<TOP  OF  FILE> 

BCRAF  *  CM70000 ,  T45 .  USER 

REQUEST/P1/*PF . 

ATTACH ( T APE3 , NLTELI NEFI LE , I D-USER  » HR- 1 > 

ATTACH( TAPE7 , BASEVIBPOP . ID-USER , HR-1 > 

ATTACH <  PROG .FASC0D2 , ID-USER , HR-1 ) 

FTN ( I -PROG , B-PROGB , L-L , T , ER  > 

EXIT(U) 

HAP. PART. 

LDSET , PRESET- I NDEF. 

LOAD.PROGB. 

SEGLOAD. 

EXECUTE. FASCOD1. 

EXIT(U) 

REUIND. OUTPUT. 

COPYCF, OUTPUT, PI. 

CATALOG ( PI , FASC0D20UT , ID-USER ) 

EOR 

FASCODE  TREE  F  ASCOD 1 - <  LA YRS »  LASER , SCANFN » PLOTT » FLTRFN » TEST ) 

FASCOD1  GLOBAL  CONSTS.FILHDR.MAIN 

FASCOD1  INCLUDE  FASCODl.BUFIN.BUFOUT 

LAYRS  TREE  LAYER**  <  PATH , HIRAC1 , HI RACQ , LBLF4 , CONTNU , ABSHRG.EHINI T 
.RADHRB.AERSOL) 

LAYER  GLOBAL  ABSORB .SCATTR.LBLF4  »L INHDR, VBNLTE-SAVE 

LAYER  INCLUDE  LAY  ER , RADFN , RDPOPS • 8K IP 

PATH  INCLUDE  PATH 

PATH  GLOBAL  PATH-SAVE 

HIRAC1  INCLUDE  H I R AC 1 . SHAPEL , SHAPES , MO I CON , RDL I N , CNVFNV , PANEL • HOLEC 

.QV.CQNT , ABSOUT .VERFN.R1PRNT ,XINT 

HIRACQ  INCLUDE  HI RACQ , SHAPEL , SHAPEG » VO I CON , RDL IN « CNVFNQ , PANELQ , 

, HOLEC , QV , CONT , ABSOUT , VERFN , RIPRNT , XINT , STRTHS 
HIRAC1  GLOBAL  FNS-SAVE 

HIRACQ  GLOBAL  F  N8Q-SAVE , XRNLTE-SAVE 

* 

LPLF4  INCLUDE  LBLF4.RDLIN4.CONVF4, HOLEC, QV, SHRINK, VOICON 

■4 

CONTNU  TREE  C0NTNH-<SLF2?A,SLF260,FRN296,FRNC02,N2C0NT> 

CONTNH  INCLUDE  CONTNH.XINT 

SLF296  INCLUDE  SLF296 

SLF260  INCLUDE  SLF260 

FRN296  INCLUDE  FRN294 

FRNC02  INCLUDE  FRNC02 

N2CONT  INCLUDE  N2C0NT 

♦ 

ABSHR8  INCLUDE  ABSMRG, ABSOUT 


* 


EHINIT 

INCLUDE 

EHINIT, EMIN.EMOUT.BBFN 

RADMRG 

INCLUDE 

RADMRG, EMIN, EMOUT ,BBFN 

SCANFN  INCLUDE  SCANFN » SHAPED  *  RDSCAN » SHRK SC , SHAPET * CONVSC , PANLSC 

,CNVREC,SINCSQ,CKPRNT 

FLTRFN 

INCLUDE 

FLTRFN » RDSCAN , CNUFL  T 

* 

TEST 

INCLUDE 

TEST 

PLOTT 

PLTFAS 

PLTFAS 

TREE 

0L09AL 

INCLUDE 

PLTFAS- (HEADER, AXEST. FPLINE) 

PLTHDR , AXISXY , YCQM , POINTS , CHEXP , TITLOC , NAME 
PLTFAS , LINT , EXPT , XNTLOG , TEMPFN , TENLOG , MNMX , L I MT 

AXEST 

EOR 

TREE 

END 

AXES-(AXISL,AXL0G,AX2) 

FASC0D1 

1  5.1  KM  SEA  LEVEL  HORZNTL 

HI-1  F4*l  CN=1  AE*0  EM*1  SC*0  FI»0  PL=0  T8*0  EQ*0 
2350.000  2400.000 
0.000 
0.000 

1 

1013.0000  288.1000  STAND.  MODEL  ATM.  6 

1.006E+23  4.254E+21  3.456E+17  3.A10E+18  9.A70E+17  2.063E+19  2.701E+24  1.007E+25 
<BOTTOM  OF  FILE> 
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APPENDIX  C  :  PROGRAM  BCDNLTE 


I 


Program  BCDNLTE  prepares  molecular  line  data  in  blocked  BCD  computer 
card  images  for  use  with  the  non-LTE  subroutines  of  FASCQD2-82.  The  pro¬ 
gram  has  been  constucted  using  program  BCDMRG  as  a  basis,  and  making  only  a 
few  modifications  so  as  to  permit  identification  of  non-LTE  transitions  by 
FASC0D2.  The  user  should  note  that  program  BCDMRG  is  completely  compat¬ 
ible  with  FASC0D2  under  circumstances  where  only  equilibrium  conditions 
(low  altitude)  are  to  be  specified.  Furthermore*  if  the  input 
parameter  INLTE  is  set  equal  to  zero,  BCDNLTE  will  function  identically 
with  BCDMRG. 


Both  BCDNLTE  and  BCDMRG  create  line  data  files  for  a  specified  spec¬ 
tral  range  for  use  by  FASC0D2.  This  data  includes  each  line's  frequency, 
intensity  and  half-width  at  296°K,  lower  state  energy  level,  and  a  molecular 
ID  between  1  and  20  which  is  used  to  distinguish  the  various  gases. 

BCDNLTE  must  also  pass  to  FASC0D2  information  to  allow  the  identification 
of  the  upper  or  lower  vibrational  states,  if  either  or  both  of  those  states 
is  nori-thermally  populated  at  high  altitude.  This  additional  information 
is  passed  through  the  last  variable  MOL,  which  in  BCDMRG  1b 
the  molecular  identifier.  In  the  non-LTE  mode,  MOL  has  the 
6-digit  integer  from  IJKLMN,  where  MN  Is  the  ordinary  molecular  ID,  1J  and 
KL  are  upper  and  lower  vibrational  state  identifiers,  respectively.  Either 
IJ  or  KL  will  be  zero  if  the  upper  or  lower  vibrational  state  doesn't  match 
one  of  the  non-equilibrium  vibrational  states,  or  if  MN  is  not  the  ID  of  a 
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II 

non -equilibrium  gas.  The  user  should  note, however ,  that  modifications  will 
also  be  required  to  FASC0D2  subroutine  RDPQPS  in  order  to  use  a  revised  set 
of  vibrational  states  in  non-LTE  calculations.  Appendix  B  provides  a 
description  of  this  routine. 

The  identification  of  the  upper  and  lower  vibrational  state  of  each 
line,  and  comparison  with  the  non-equilibrium  states  of  the  model  of  Degges 
and  Smith  is  accomplished  in  subroutine  VIBQU.  This  routine  is  called  orice 
for  each  line.  Standard  vibrational  ntato  ID's  are  compared  against  the 
specific  non-equilibrium  state  ID's  which  have  been  stored  in  array  DATA 
statements.  If  a  match  occurs,  the  corresponding  number  (presently  from 
1-10)  is  added  to  the  proper  digits  of  MOL  (forming  M0L2  from  M0L1  in 
VIBQU).  The  number  of  non-equilibrium  vibrational  states  is  currently 
taken  to  be  B  for  water  vapor,  10  for  carbon  dioxide,  2  for  ozone,  and  3 
for  nitric  oxide.  In  each  case,  the  ground  state  is  included  first. 

Updating  and  modification  of  the  selected  set  of  vibrational  states  for 
each  gas  can  be  done  simply  by  modifying  the  two  DATA  statements  within 
VIBQU  which  initialize  NUMH20 ,  NUMC02,  NUM03 ,  NUMNO ,  and  HVQH20 ,  HVQC02, 

HVQ03,  HVQNO. 
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C.l  SAMPLE  FILE  FOR  THE  SUBMISSION  OF  BCDNLTE 


<T0P  OF  FILE> 

SCRAF»T300*HT2.  # m 

REQUEST /BCDOUT /TAPE3/*PF . 

ATT  ACH  ( PROG ,  BCDNLTE ,  I  D*USER ) 

FTN » I *PR0G . B»PROGB » T  *  £R . 

VSN(TAPE1*CC0157) 

LABEL(TAPE1*R*L*CC0157) 

VSN ( T APE2*CC0 171) 

LABEL(TAPE2»RiLsCC0171) 

C0PYCF*TAPE1*DUM*22. 

COP YCF * T APE2 * DUN  *  22 . 

PROGB. 

EXIT(U) 

CATALOG ( TAPE3 , NLTELINEFILE* ID*USER ) 

EXIT(U) 

REWIND* OUTPUT. 

COPYCF, OUTPUT* BCDOUT. 

CATALOG (BCDOUT  *BCDNLTEOUT»ID*US£R) 

RETURN(  TAPED 
RETURN (TAPE2) 

EXIT  <U> 

EOR 

RUN  OF  BCDNLTE  FOR  USE  IN  NLTE  CALCULATIONS 
2300.000  2450.000  I NLTE*  1 

niniiiuuiiniiii  tape* 

limiiuiumum  trace 

<BOTTOM  OF  FILE> 


86 


USER 


nnnnn  nonnonnnon 


C.2  PROGRAM  BCDNLTE :  LISTING  OF  SUBROUTINE  VIBQU 


SUBROUTINE  VIBQU(LTEHP) 


THIS  SUBROUTINE  COMPARES  THE  TUO  VIBRATIONAL  LEVEL 
ID'S  OF  SELECTED  MOLECULES  AGAINST  STORED  VIB  STATE 
ID'S.  MOL  IS  ADJUSTED  TO  INCLUDE  ADDITIONAL  IDENTIFIERS. 

IJ - UPPER  STATE  ID 

KL - LOWER  STATE  ID 

MN  -  ORDINARY  MOLECULAR  ID 


COMMON  TI7(7)tTI10(10)«R0*R(320) 

DIMENSION  HVQH20 <  8 ) » HVQC02 < 1 0 ) » H VQ03  <  2 ) » H VQNO ( 3 ) 
EQUIVALENCE  <TI7<7> ,M0L2> 

EQUIVALENCE  UI10<5>,VIB1MTI10<A>,VIB2MTI10<10>,M0L1> 
DATA  NUMH20 « NUMC02  *  NUM03 • NUMNQ/8  *10*2*3  / 

DATA  (HVQH20<IM*1,8>/ 

1  AH  0  0  0, 

2  AH  0  1  0, 

3  AH  0  2  0, 

4  AH  1  0  0* 

5  AH  0  0  1. 

A  AH  0  3  0, 

7  AH  1  1  0, 

8  AH  0  1  1  / 

DATA  (HVQC02(I >.1*1,10)/ 

1  10H  0  0  0  0  1, 

2  10H  0  1  1  0  it 

3  10H  10002i 

4  10H  0  2  2  0  1, 

5  10H  11102* 

A  10H  0  3  3  0  1* 

7  10H  1  0  0  0  1, 

8  10H  0  0  0  1  1« 

9  10H  10012* 

C  10H  1  0  0  i  1  / 

DATA  < HVQ03 <I)»I=1*2>/ 

1  AH  0  0  0. 

2  AH  0  0  1  / 

DATA  (HVQNO(I) i 1*1 »3)/lH0,lHl,lH2  / 


M0L1  IS  THE  ORIGINAL  MOLECULAR  ID 
M0L2  IS  THE  NEW  TAGGED  MOL  ID 

H0L2=M0Li 
C 

IF(NOLl.NE.l)  GO  TO  200 
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C  WATER  VAPOR  LINE 
C 

REWIND  LTEHP 

WRITE(LTEHP»900)  <TI10<J),>5,9) 

C 

REU1ND  LTEHP 

READ (LTEHP, 901)  VIB1,VIB2 

900  FORMAT <5A7> 

901  F0RMAT(21X,A6,1X«A6) 

NUPP=0 

DO  150  NVQ*1,NUMH20 
HTEST=HVQH20(NVQ) 

IF  (VIB1.NE.HTEST)  00  TO  ISO 
NUPP-NVQ 
GO  TO  160 
150  CONTINUE 
160  NL0W“0 

DO  170  NVQ=1»NUMH2Q 
HTEST ■HVQH20 ( NVQ ) 

IF  ( VI B2.NE. HTEST)  00  TO  170 
NLOWsNVQ 
00  TO  180 
170  CONTINUE 

180  H0L2  «  H0L1  +  100*NLOU  +  10000*NUPP 
00  TO  500 
200  CONTINUE 

IF<M0L1.NE.2)  00  TO  300 
C 

C  C02  LINE 

C 

REWIND  LTEHP 

WRITE  (LTEHP,  900)  <T110(  J)  ,>5,9> 
REWIND  LTEHP 

READ (LTEHP, 902)  VIB1,VIB2 

902  F0RMAT(2X,A10,5X,A10> 

NUPP«0 

DO  250  NVQ*1,NUHC02 
HTEST “HVQC02 ( NVQ ) 

IF  (VIBl.NE. HTEST)  GO  TO  250 
NUPP*NVQ 
00  TO  260 
250  CONTINUE 
260  NIOW-O 

DO  270  NVQ=1,NUHC02 
HTEST =HVQC02 ( NVQ  > 

IF  (VIB2.NE. HTEST)  00  TO  270 
NL0W=NVQ 
80  TO  280 
270  CONTINUE 


280  MDL2  =  M0L1  +  100*NL0U  +  10000+NDPP 
GO  TO  500 
300  CONTINUE 
C 

IF(H0L1.NE.3)  00  TO  400 
C 

C  OZONE  LINE 
C 

REWIND  LTENP 

WRI TE( LTENP » 900)  <TI10<J),J*5,9) 
REWIND  LTENP 

READ (LTENP? 903)  VIB1?VI82 

903  FORMAT < 21 X,A4,iX,A4> 

NUPP*0 

DO  350  NVQ»1?NUN03 
HTEST*HVQ03<NVQ) 

IF  (VIB1.NE.HTEST)  00  TO  350 
NUPP-NUQ 
GO  TO  340 
350  CONTINUE 
340  NLOU-O 

DO  370  NUQ-1.NUH03 
HTEST “HVQ03 ( NVQ ) 

IF  (VIB2.NE. HTEST)  00  TO  370 
NLOU-NUQ 
00  TO  380 
370  CONTINUE 

380  H0L2  -  NOLI  +  100*NL0W  4  10000*NUPP 
GO  TO  500 
C 

400  CONTINUE 

IF(MOLl.NE.B)  GO  TO  500 
C 

C  NITRIC  OXIDE  LINE 
C 

REWIND  LTENP 

URITE(LTEMP,90Q>  <TI10<J),.:*5,9) 
REWIND  LTEHP 

READ (LTENP? 904)  VIB1.YIB2 

904  FORMAT (7X?A1?7X?A1) 

NUPP*0 

DO  450  NVQ*1 ?NUMNO 
HTESTsHVQNO<NVQ) 

IF  (VIB1.NE. HTEST)  GO  TO  450 
NUPPsNVQ 
GO  TO  440 
450  CONTINUE 
440  NL0W=0 

DO  470  NVG=1 iNUMNO 
HTEST=HVQNO(NMQ) 


IF  (VIB2.NE.HTE3T)  QO  TO  470 
NLOU»MVQ 
00  TO  480 
470  CONTINUE 

480  H0L2  «  M0L1  4  100*NL0U  4  10000*NUPP 
300  CONTINUE 
C 

DO  550  J*l,6 
TI7<J)-TI10(J) 

550  CONTINUE 
RETURN 
END 


so 


